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ABSTRACT 


This IS Volume 2 of three volumes which report the results of a strapdown calibration 
and alignment study performed by the Univac Federal Systems Division for the Guidance 
Laboratory of NASA/ERC. 

This study develops techniques to accomplish laboratory calibration and alignment of 
a strapdown inertial sensing unit (ISU) being configured by NASA/ERC. Calibration is 
accomplished by measuring specific input environments and using the relationship of 
known kinematic input to sensor outputs, to determine the constants of the sensor 
models. The environments used consist of inputs from the earth angular rate, the 
normal reaction force of gravity, and the angular rotation imposed by a test fixture 
in some cases. Techmques are also developed to accomplish alignment by three 
methods. First, Mrror Alignment employs autocollimators to measure the earth 
orientation of the normals to two mirrors mounted on the ISU. Second, Level Align- 
ment uses an autocollimator to measure the azimuth of the normal to one ISU mirror 
and accelerometer measurements to determine the orientation of local vertical with 
respect to the body axes. Third, Gyrocompass Alignment determines earth alignment 
of the ISU by gyro and accelerometer measurement of the earth rate and gravity 
normal force vectors. 


The three volumes of this study are composed as follows. 


• Volume 1 — Development Document. This volume contains the detailed 
development of the calibration and alignment techniques. The develop- 
ment IS presented as a rigorous systems engineering task and a step- 
by-step development of specific solutions is presented. 

• Volume 2 — Procedural and Parametric Trade-off Analyses Document. 
This volume contains the detailed trade-off studies supportii^ the 
developments given in Volume 1. 

• Volume 3 - Laboratory Procedures Manual. In Volume 3 the imple- 
mentation of the selected procedures is presented. The laboratory 
procedures are presented by use of both detailed step-by-step check 
sheets and schematic representations of the laboratory depicting 
the entire process at each major step in the procedure. The equa- 
tions to be programmed in the implementation of the procedures 
are contained in this volume. 
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GLOSSARY 


As an aid to understanding the symbolism, we present the following rules of notation. 


• Wherever possible symbols will be used which suggest the name of the 
parameter involved. 

o Lower case subscripts are used almost exclusively for indexing over several 
items of the same kind. Examples are the indexes used to identify the three 
gyros, the three accelerometers, the two pulse trains of each accelerometer, 
the two clock scale factors, etc. 

• Lowercase superscripts are used to index over different positions. 

• Uppercase superscripts and subscripts will be used to distinguish between 
parameters of the same kind. For example, T is used to identify a 

RE 

transformation matrix. Lettered superscripts such as BE in identify 
the particular transformation. 

• An underline will identify a vector. 

• Unit vectors are used to identify lines m space such as instrument axes 
and the axes of all frames of reference. 

• Components of any vector along with any axis is indicated by a dot product 
of that vector with the unit vector along the axis of interest. 

• The Greek sigma (E) will be used for summations. Where the limits of 
summation are clear from the context, they will not be indicated with the 
symbol. 

• The Greek A is always used to indicate a difference. 

• S 0 and C <6 are sometimes used to identify the sine and cosine of the angle 0. 

• A triple line symbol (=) will be used for definitions. 

• A superior " " denotes a prior estimate of the quantity. 

• A superior" ^ "denotes an estimate of the quantity from the estimation routine. 
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(Aj* B^) 


Ax 




Applied acceleration vector. 

Elements of 

Unit vector directed along the input axis of the ith accelerometer 
i = 1, 2, 3. 

A vector determined by the Alignment Parameter Evaluation 
Procedure and input to the Estimation Houtine. 

Unit vector directed along the ith Body Axis i = 1, 2, 3. 


Bp Bq, Bg Gyro unbalance coefficients. 

Compliance Coefficients. 


Counters 


D 

D 

D 


0 

1 

2 


The six frequency counters used as data collection devices 
during calibration. 

Accelerometer bias. 

Accelerometer scale factor. 

Accelerometer second order coefficient . 


D 


3 


Accelerometer third order coefficient. 


E Unit vector directed East (E^). 

E . Unit vector directed along the ith Earth Axis . 


Eq Quantization error. 

L , fo Frequencies of accelerometer strings 1 and 2, in zero 

crossings per second. 


F 

-1 


A triad of orthogonal unit vectors attached to the base of the 
table. 


Gj Unit vector directed alor^ the ith input axis of the gyro. 

Elements of 


g The vector directed up that represents the nornoal force to 

counteract gravity in a static orientation. Corresponding to 
popular convention, this is referred to as the "gravity vector". 

1/ O Input/ Output. 


x 



lEU 

ISU 

J 

K 


m 

M 


Ml 

N 

NpN2 

"g 

En| 

o 



Triad of orthogonal unit vectors attached to the inner axis of 
test table. 

Interface Electronics Unit - system interface device for the 
laboratory computer. 

Inertial Sensing Unit. 

Gyro angular rate coefficient. 

Number of samples of accelerometer and gyro data taken in 
Alignment. 

Position index used in calibration (superscript) . 

Matrix generated by Alignment Parameter Evaluation and used 
by Alignment Estimation Routine. 

Unit normal to ith mirror. 

Unit vector directed North (Eg). 

Count of output pulses from strings 1 and 2 of accelerometer. 
Instrument noise in accelerometer. 

Instrument noise in gyro. 

Count of output pulses from strings 1 and 2 of accelerometer. 

Count of timing pulses from master oscillator to frequency 
counters. 

Count of timing pulses from master oscillator to lEU. 

Unit vector directed along the output axis of gyro. 

Triad of orthogonal unit vectors attached to the outer axis 
of the table. 

Unit vector in the direction of the projection of M- in the plane 
formed by E and N. ^ 

Defined on Chart 4-12 of the Development Document. 

Defined on Chart 4-4 of the Development Document. 

The transformation from accelerometer input axes to body axes." 
The transformation from gyro input axes to body axes. 


XI 



Qjp Qjs 


r 

R 




Resolver 




t 

T 


rj^Bl 

,pBRm 


U 

V 

W 

X-Y 


Z 


(Erf,) 

(S6)j 


Gyro dynamic coupling coefficients. 

Position vector. 

Gyro bias. 

Triad of orthogonal unit vectors attached to rotary axis of 
table. 

Angular resolvers on each axis of the test table. 

Unit vector directed along the ith gyro spin axis. 

Scale factor associated with pulsed output from test table rotary 
axis. 

Scale factor associated with timing pulses accumulated by the 
frequency counters. 

Scale factor associated with timing pulses to the lEU. 

Time. 

In alignment, the determined alignment matrix to transform 

SE 

from body to earth axes. T is equivalent to T . 

Transform from ISU body axes to inner axis frame. 

Transform from ISU Body Frame Axes to Rotary Axis Frame in 
the mth orientation. 

Triad of orthogonal unit vectors attached to the trunnion axis 
of the test table. 

Unit vector directed up (Ej^). 

Velocity vector. 

E 

Unit vector directed along w . 

Dual input on frequency counter that will difference two pulse 
trains for comparison with a third input (Z). 

Input on frequency counter for pulse train. 

The azimuth angle of the normal to the ith mirror. 

Pulsed output from the jth string of the ith accelerometer. 

Pulsed output of the ith gyro. 

Gyro scale factor. 
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"c 

The clock quantization error. 

€t 

In instantaneous alignment estimation techniques, this symbol 
represents the length of time after completion of the last 
measurement to the time at which the prediction is made. 


The zenith angle of the normal to the ith mirror. 

X 

Local colatitude. 

a 

g 

The estimated rms error in the magnitude of g. 

^0 

The estimated rms error in the direction of "up". 

01 

Angular displacement about the trunnion axis of the test 
table for calibration position m. 

^2 

Angular displacement about the rotary axis of the test table for 
calibration position m. 

^3 

Angular displacement about the outer axis of the test table 
for calibration position m. 

< 

Angular displacement about the inner axis of the test table 
for calibration position m. 

0n(t) 

Covariance function of accelerometer noise used in Alignment 
Parameter Evaluation. 

0 (t) 

Covariance function of translational acceleration noise used 
in Alignment Parameter Evaluation. 

0{t) 

Covariance function of rotational noise used in Alignment 
Parameter Evaluation. 

w 

Angular velocity vector. 

T 

CO 

Angular velocity of the test table rotary axis. 

E 

w 

Earth rotation vector. 

E 

CO cos X 

Component of earth rotation vector alot^ the vertical. 

E 

00 sin X 

Component of earth rotation vector along north. 
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SECTION 1 
INTRODUCTION 


This document, in conjunction with two other volumes, describes the achievements of a 
six month study conducted for the: 


Guidance Laboratory 
Electronics Research Center 
National Aeronautics and Space Administration 
Cambridge, IVIassachu setts 


by the: 


Aerospace Systems Analysis Department 
Univac Federal Systems Division 
Saint Paul, Minnesota 
A Division of Sperry Rand Corporation 

The purpose of the study is to develop techniques and outline procedures for the laboratory 
calibration and alignment of a strapdown inertial sensing unit. The Development Docu- 
ment, Volume 1, presents a detailed analysis of the calibration and alignment problem 
and develops a specific solution. This document, Volume 2, is a set of addenda which 
serve to justify the conclusions reached in the development of specific calibration and 
alignment techniques in Volume 1. Reference is made to the Development Document 
throughout the presentation of the procedural and parametric trade-off analyses. The 
Laboratory Procedures Manual, Volume 3, describes the procedures for an operational 
implementation of the solutions obtained in Volume 1. A statement of the study objectives 
is contained in Section 1 of the Development Document. 

The trade-off analyses described herein assume a laboratory facility which includes 
a test table, a computer and interface unit, frequency counters, autocollimators, three 
gyroscopes and three accelerometers. This laboratory facility is described in detail 
in Section 3 of the Development Document. The instruments are described in detail in 
Section 2. 2 of the Development Document. This laboratory defimtion was not considered 
an absolute constraint on the trade-off studies, however, for the effect of variations in 
laboratory equipment and instrument characteristics have been carefully considered. The 
intent of the trade-off analyses has been to provide sufficient data on the relationships 
between important trade-off parameters to assist the laboratory test program at ERC. 

This program is directed toward the larger problem of laboratory testing a strapdown 
inertial sensing unit for the purpose of evaluating many advanced guidance and navigation 
concepts. 
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It is beyond the scope of the study, however, to consider all of the parameters which 
would be important in the ISU application to any specific mission. Some of the interesting 
parameters could not reasonably be considered without a great deal of presently unavail- 
able supporting data. For example, parameters such as instrument reliability, instru- 
ment cost, instrument location relative to the body and to each other, etc. 

The parameters that have been considered in the trade-off analyses were agreed to by 
NASA/ERC and do provide a wealth of information on the variables one would like to 
control in designing an ISU for space application. These trade-off parameters have been 
organized into the trade-off analyses reported in this document. 

Figure 1-1 (or 1-2) illustrates the dependence of calibration (or alignment) time and 
accuracy upon the selected parameters. Note that all of them are related to either . 
accuracy or time for both calibration and alignment. Thus if either accuracy or time is 
constrained one would have a solvable optimization problem. Given a set of ISU accuracy 
requirements, one could select the instrument characteristics (in terms of calibration 
coefficient stability, internal noise, and readout quantization), environmental noise con- 
straints required, alignment time, and calibration time which would satisfy the require- 
ments in a "best" way. As used here, "best" implies minimization of instrument design 
and production complexity, calibration time, and alignment time. 

The specific calibration and alignment procedural and parametric studies described in 
this document are listed by general category in Chart 1-1. The sections in which the 
studies are covered are listed in the margin of this chart. 

The calibration trade-off analyses are the subject of Section 2. The alignment trade-off 

analyses are described in Sections 3, 4, 5, and 6. Section 3 serves as the introduction 

to alignment trade-offs. In Section 4, the general error equations relating errors in the 
IS 

estimate of g and w to a basic measure of alignment error are developed for all four 

orientations and for both level and gyrocompass alignment. The trade-offs leading to 

E 

the selection of the g and w estimation techniques are the subject of Section 5, The 
expected alignment errors for these techniques as a function of instrument and environ- 
ment noise, instrument readout quantization, sample time and number of samples, 
estimation iteration, and computer word length are also developed in Section 5 from a 
computer simulation. Section 6 develops alignment accuracy as a function of calibration 
accuracy. Both worst-case and one sigma errors are treated. 

As a matter of easy reference, we list below a cross reference between the trade-off 
studies called for in the Statement of Work, and those covered in this document: 
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Figure 1-1. Calibration Trade-Off Parameters 

















Figure 1-2. Alignment Trade- Off Parameters 

















CHART 1-1 


Calibration, and Alignment Trade-off Studies 


Calibration Trade-offs 


Procedural Trade-offs 

Calibration Time vs Calibration Accuracy 

A Because of instrument quantization 
A Because of instrument and environment noise 


Alignment Trade-offs 


3 Alignment Time vs Alignment Accuracy 
A Because of Instrument quantization 

• Alignment Matrix Accuracy vs Precision of the 
Estimation of Body Axes Components of Gravity 
and Earth Rate 


• Estimation of Body Axes Components of Gravity 
and Earth Rate vs Estimation Technique 

A Sampling rate considerations 
A Total sampling time considerations 
A Word length 
A Algorithm 

• Alignment Accuracy vs Calibration Accuracy 


Section 

Number 


L-5 








1. Alignment Time versus Alignment Accuracy 
Because of Worst-Case Quantization — Section 3. 3 
Because of Simulated Noise and Quantization — Section 5 

2. Alignment Time versus Sensor Quantization 
Worst-Case — Section 3. 3 

Simulated Quantization - Section 5 

3. Alignment Time versus Iterative Scheme — Section 5 

4. Alignment Accuracy versus Word Length — Section 5 

5. Alignment Accuracy versus Data Sampling Rate — Section 5 

6. Alignment Accuracy versus Calibration Accuracy — Section 6 

7. Calibration Time versus Calibration Accuracy - Section 2.2 

8. Calibration Time versus Calibration Procedure — Section 2. 1. 


1-6 



SECTION 2 



SECTION 2 

CALIBRATION TRADE-OFF STUDIES 


Not all of the calibration trade-off analyses can be expressed analytically. Certain of the 
trade-offs are procedural in nature and deal with variables which cannot be measured 
quantitatively. Such variables were encountered as a choice to be made between a small 
number of alternatives before one could proceed in the calibration technique develop- 
ment. These procedural trade-offs are discussed in detail in the first subsection. 

Calibration time is related to accuracy through the effects of instrument internal noise, 
instrument readout quantization, and environmental noise. The parametric analyses of 
calibration time and accuracy as a function of these parameters is the subject of the 
second subsection. Plots of calibration coefficient accuracy versus data collection time, 
as a function of worst- case quantization and statistical noise, are included at the end of 
that section. These plots are reproduced in Section 11-6 of the Laboratory Procedures 
Manual. 

2. 1 PROCEDURAL TRADE-OFFS 

Throughout the developments of the calibration techmques in Volume 1, procedures are 
chosen where alternatives could have been taken. In most instances explanations are 
given as to why the choices taken are better than the alternatives. In order to preserve 
the smoothness of presentation, the explanations are always rather brief in the Develop- 
ment Document. In this section of the trade-off document we present more detailed 
explanations for the more important decisions made in the Development Document. The 
presentation takes the form of a listing of explanations, where the order more or less 
indicates the relative importance. 

2. 1. 1 Isolation of Calibration Constants 

In the introduction to Section 4. 2 of the Development Document it is mentioned that the 
calibration of the n constants in any instrument equation can be accomplished by the 
simultaneous solution of n calibration equations, where each equation corresponds to 
the input/ output relationship for a particular environment. This technique involves, in 
general, the inversion of a n x n matrix, which would be a very cumbersome approach 
to calibration. 
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Our approach has been to choose environments such that a large number of the unknown 
parameters are not sensed by the instruments. This approach has several advantages 
over the more general inversion technique. The most important advantage is that a use 
of environments where the general calibration equations are considerably reduced allows 
for a satisfaction of the precision requirements on only those constants which are sensed. 
In environments where all terms are sensed, it is impossible to reduce the equations, 
and therefore it is impossible to control the precision on all constants independently. 
Another advantage of the use of reduced equations is that large matrices do not have to 
be inverted. Operational matrix inversion routines are always approximations when the 
matrices are large, as they would be for a general use of our equations. In the chosen 
techniques described in Section 4 of the Development Document, we never use more than 
two equations in the solution of any constants. Therefore, the inversions always involve 
2x2 matrices, which can be accomplished with no approximations. 

2. 1. 2 Use of Frequency Counters for Data Collection 

In Section 4. 3 of the Development Document it is mentioned that frequency counters (see 
Section 3. 3) are to be used for all instrument data collections. For the defined ERC 
facility there is only one alternative to the use of frequency counters, that being a use of 
the laboratory computer interface. The frequency counters have one major advantage 
over the computer interface, that being the ability to cycle off of the instrument data 
train (that is, the ability to detect the leading edge of the instrument pulses), as opposed 
to cycling off of a clock. In reading the accelerometer outputs this is a distinct advan- 
tage, for the quantization error would be merely the uncertainty in reading the leading 
edge of the pulses, which was assumed to be one-tenth of a pulse for each of the 
accelerometer pulse trains. The computer interface, on the other hand, can only be 
cycled off of clock time, and the worst- case quantization error would be one pulse per 
pulse train. In reading the gyro outputs, the advantage of using the frequency counters 
is not so great as for the accelerometers. The gyro torquing is driven by a clock, re- 
sulting in the leading edges of the pulses not occurring simultaneously with increments 
of angular rotation. Such pulse trains have an implicit worst-case error of one quantum. 
Reading the leading edge does not reduce the error by any amount. A sampling of the 
gyro pulse trains asynchronously in time with the computer could result in a worst-case 
error of as much as the two quanta, so there is a small advantage in using the frequency 
counters. The principal reason for using the frequency counters for collecting the gyro 
data is not this small advantage. The principal reason is dictated by the necessity for 
keeping the computer out of the data collection process so that computer malfunctions 
would not disrupt the calibration. 
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2.1.3 Calibration of Gyro and Accelerometer Acceleration Sensitive Terms in Static 
Positions 


Recall that we mention m the Development Document that the reason for not rotating the 
table during accelerometer calibration is so that no rotation-induced accelerations 
would be introduced to the accelerometers. All such induced acceleration environments 
would have to be measured independently for the purpose of calibration. Those rotation- 
induced accelerations would be functions of the distance from the axis of rotation, and it 
would be difficult to ascertain the radius which locates the sensitive point in the accel- 
erometers. Additionally, the rotationaUy- induced accelerations would not be sufficiently 
large to be considered useful as inputs. There is, therefore, no reason to rotate the 
table during accelerometer calibration. 

The reason for not rotating the gyros during the calibration of acceleration-sensitive 

O 1 

terms is to minimize the influence of the scale factor and (Q )~ matrix imprecisions on 
the calibration of unbalance and compliance coefficients. 

2.1.4 Use of Test Table Rather than Autocollimators for Measurement of Environment 
In Section 4. 2 of the Development Document the body-axis-coraponents of g, o) > and 

tji — — 

(0 are expressed in terms of test table parameters. In Section 3. 2. 1. 3 of the same 

F' 

document it is stated that the transform of g and w from earth axes to body axes can be 
accomplished with the assistance of the autocollimators as well as the table resolvers. 

As a matter of fact the autocollimators might transform the vectors more accurately, 
for the resultant transformation would involve a smaller number of error sources. Un- 
fortunately the transformation to body- axes- components must be accomplished for many 
different orientations of the ISU. To use the autocollinoators for all positions would re- 
quire a great deal of time. For each new position, the optical devices would have to be 
moved to a different location in the laboratory and resurveyed. Either that, or a large 
number of autocollimators would have to be purchased, two for each nominal position. 

A small compromise on precision must therefore be accepted, and the test table used 
for all transformations. 

2.1.5 Whole Turn Data Taken During Calibration of Gyro Scale Factors and 
Misalignments 

It is possible to calibrate the gyros in the presence of table motion without taking data 
from whole turns. The reason for using whole turn data only is that it is the simplest 
and most accurate technique. Taking data from fractions of whole turns would introduce 
the following undesirables: 
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• Fractions of whole turns of the test table cannot be measured as accurately as 
whole turns. That is, 360° can be measured more accurately, with a resolver 
than, say, 20° or any other fraction of a whole turn. 

• The transient terms (integrals of sine and cosine of 0^ or 202) m Chart 4-8 of 
the Development Document will always be evidenced in fractions of whole turns. 
The evaluation of those integrals would have to be accomplished digitally, which 
would introduce errors into the calibration. 

• The whole turn equations allow for the nulling of linear acceleration inputs to 
the gyros (that is unbalance terms). Those terms can be nulled over whole 
turn integrations only. 


It must be noted that the last two advantages would not be available if it were not for the 
fact that the table speed can be regulated sufficiently close to a constant. 

2. 1. 6 Use of Maximum Speed of Table During Scale Factor and Misalignment Calibration 


It was mentioned in Section 4. 2. 1 of the Development Document that the table will always 
be rotated at the maximum allowable speed (below the saturation of the gyros) during the 
data collection from the first six positions. The reasons are very simply stated: 


• For a given precision requirement on the constants, the higher the speed the 
less the amount of time required for data collection. This is evidenced in the 
calibration equations found in Section 4. 3. 2 of the Development Document. 

• During calibration of w- sensitive coefficients it is advisable to make the 
angular velocity terms predominate. This is accomplished by using the highest 
possible speed. 


Because the gyro scale factor is unknown at the beginning of calibration, the saturation 
level is not precisely known. Therefore, the speed will have to be regulated, experi- 
mentally, to be close to but be less than gyro saturation which introduces at least one ' 
pulse sign change over a finite period of time (say 10 seconds). 


2. 2 CALIBEATION TIME VERSUS CALIBRATION ACCURACY 


In the development of the calibration equations (Section 4, Development Document), the 
effects of internal noise, environmental noise, and readout quantization on instrument 
output were neglected. Thus, the gyro output 


— S 6r 


Eq 
A 4> 


An 

A 4> 


Ai« ^ 
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was approximated by 


P 


G 


S6k 
k ^ 


so that we define 


AP^ = — [eq + An + f^^Aco.Gdtl 
A^ 0 


( 1 ) 


In the same way, the accelerometer output (where internal noise is assumed negligible) 

= (Ng - N^) - Ell - D^J^^^^Aa.Adt 


was approximated by 


P^ = (Ng - N^) 


so that we define 

AP^ = Eq + D^J^^^Aa-Adt (2) 

where An is the internal noise of the gyro, 

Aw and Aa are the environmental noise, 

Eq is the quantization error. 

These approximations produce errors in the calibration parameters which are functions 
of data collection time. In this subsection then, calibration accuracy will be studied as 
a function of calibration time. 

The general calibration parameter error equation is developed in Section 2. 2. 1. The 
specific error equations for quantization and noise are developed in Sections 2, 2, 2 and 
2. 2. 3, respectively. Finally, the calibration parameter errors are plotted as a func- 
tion of data collection time. At, in Section 2. 2. 4. 

It should be noted that while the errors plotted due to noise are statistical (standard 
deviation), the errors plotted due to quantization are worst case. The errors due to 
noise and quantization are treated as independent in the development which follows. 
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Thus, noise-free instruments and environment are assumed in the study of readout 
quantization, and conversely. 

2. 2. 1 General Error Equation 

The equations for calibration parameters, as a function of gyro or accelerometer data 
collected, are derived in Section 4. 3, 2 of the Development Document. All of these 
equations are linear and have the form 


P 

y = A ~ + B (3) 

At 


where y is the parameter to be calibrated 
A IS .some known constant 
B is a function of previously computed quantities 
At is the data collection time interval 
P IS the instrument readout summed over At. 

Since equation (3) is linear, the error In y due to an error in P can be expressed by 

A 

Ay = AP (4) 

At 

In the two subsections which follow, the AP due to noise, and the AP due to quantization 
of the instrument readout are discussed. 

2.2.2 Quantization Errors 

The magnitude of the error in the instrument readout due to quantization of the readout 
varies with the particular instants in time at which readout is commenced and ended. 

A worst- case error has been assumed to be 

AP*^ = 1 pulse 


for the gyro, and 


AP'^ = 0. 1 pulse 


for the accelerometer. 
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Equation (4) can be rewritten in the form, with Ay = 

jS-At = A(AP) 1 = 1, 2, 3, 4 (5) 

The plot of the worst-case ySj against At is a family of h 3 q)erbolas. Chart 2-1 which 
follows, identifies the quantities ^8^^ for the gyro. Chart 2-2 identifies the quantities 
for the accelerometer. 

Figures 2-1 through 2-4 show plots of the jS^^'s for worst-case quantization. 

2.2.3 Noise Errors 

The analysis of noise- induced errors is based on equations (1) and (2) in the introduction 
to 2. 2 after dropping the quantization terms. In equation (1) the integrand Aw> G includes 
a component of earth angular rate modulated by an angular displacement noise, A9, in 
addition to the angular velocity from environment disturbances, which is sensed directly 
by the gyro. The development in Appendix A shows that the error input to the gyro is: 

Aw* G = (o)^ sin e^) AS + 

E E 

where co is the earth rate (magnitude); 9 = 0^ + AS is the angle between w and _N 

(colatitude); and w® is the environment angular velocity noise, the time derivative of A0. 

The environment angular displacement noise AS (and hence also its derivative, w®) is 

assumed zero about a vertical axis and isotropic in the horizontal plane. Therefore: 

A«* G = w®(A9/^/T ) + w®, for gyro horizontal 
E . r- 

= w A0 / V 2 , for gyro vertical 
where 9^^ has been assumed 45°. 

Similarly, for the accelerometer, the integrand in equation (2) Aa* A includes a component 
of gravity modulated by the angular displacement noise, A0; in addition to a®, the accelera' 
tion from environment disturbances. As developed in Appendix A, the error input to the 
accelerometer is: 


Aa* A = a® + (g sin 0^) A9 
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CHART 2-1 


GYRO COEFFICIENT ERRORS FROM WORST CASE QUANTIZATION 


1) P-.At = — (deg) 
^ 300 


2) ^„At = 


= 


(1/A^) At 


(deg/hr) 


= AR, the gyro bias term 
= ABa, the gyro unbalance term 
= ACa^, the gyro compliance term 


c 7(1/A^>) At 


(dimensionless) 


A(l/Aa>) 

= ^ the gyro scale factor term 

(1/A«) 

= A(G.* Bj), the gyro alignment term 
= AQct) the gyro nonlinearity term 


5) = 43200 (deg/hr) 

6) A* =12 (sec/pulse) 

7) AP = 1 (pulse) 
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CHART 2-2 


ACCELEROMETER COEFFICIENT ERRORS FROM WORST CASE QUANTIZATION 


1) ^gAt=- 


2) jS.At = (g-sec) 

^ 1270 


(gD. )At 


(dimensionless) 


the accelerometer scale factor term 


= A(Aj^* the accelerometer alignment term 

= ADga, the accelerometer second order terra 

= ADgR^, the accelerometer third order term 
AP 

^4= (g’s) 

^ D^At 

= ADg, the accelerometer bias term 


a = g 


4) = 254 (pulse/sec/g) 

5) AP = 0. 2 (pulse) 
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where is the magnitude of the environment acceleration noise, g is the (constant) 
acceleration of gravity; and 0 = 6^ + A0 is the angle between A and the local vertical. 
Environment acceleration noise is assumed isotropic. Since 0^ is different in the 
horizontal and vertical orientations, we have: 

Aa* A = a® + g AS , for accelerometer horizontal 
e 

= a , for accelerometer vertical 


Thus for both gyro and accelerometer the effect of environmental noise is strongly de- 
pendent upon orientation of the instrument. 


Substituting equations (6) into (1), and equations (7) into (2), and dropping the quantization 
terms, we get expressions for AP needed to implement equation (5) for the noise case. 
Since noise errors are indeterminate, a statistical approach is used, giving an rms value 
for the coefficient errors, 0^. The resulting equations are developed in Appendix A, and 
are listed in Cliarts 2-3 and 2-4. The functional relationships between noise- induced 
coefficient errors and calibration time are plotted in normalized form in Figures 2-1 
through 2-4. The rest of this subsection outlines the approach and assumptions made: 

o The variance (a ) of AP is calculated from the power spectral density of AP, 
which in turn is found from the power spectral densities of the input noises. 

• Gyro internal noise* is introduced as an equivalent noise input. 

• Accelerometer internal noise is negligible. 

• Sensor dynamics are effectively neglected. That is, the sensor transfer function, 
T{s), is assumed flat out to a frequency where the noise falls off drastically, 

• Intersample time is neglected; i. e. , we assume continuous data instead of 
sampled data. This means replacing the summation of pulses with an integration. 

4 

• All noise sources are assumed stationary, independent, random processes of 
zero mean. In addition, angular displacement noise is assumed isotropic in the 
horizontal plane, and zero in the vertical. Acceleration noise is assumed 
isotropic. 

• Environment noise is taken from Section 3. 2. 2 of the Development Document. ** 
The acceleration environment noise spectrum, Pa(f), is given in Figure 3-3 of 
that subsection, and the angular displacement noise spectrum (tilt), Pe(f), is 
given in Figure 3-4. The angular velocity noise spectrum, Pj^(f), is obtained 
from Figure 3-4; by multiplication by ( 27 rf)^. 


*GG334 Gas Bearing Gyro, Technical Description ASD-3, Honeywell Aeronautical 
Division, Minneapolis, 28 November 1966. 

**Also from H. Weinstock, Limitations on Inertial Sensor Testing Produced by Test 
Platform Vibrations, NASA TN D-3683, Washington, D. C. , November 1966. 
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CHART 2-3 


GYRO COEFFICIENT ERRORS FROM STATISTICAL NOISE 




1/2 


[j ' 

AtL’*0 


Pp(f) df 


1/2 


2) 

where 

1) 


2 ) 


^^2^rms ' ^ ^^l^rms 

oy 


AP 

= 7 , " T (deg/hr) 


(l/A<S)At 

= AR, the gyro bias term 

= ABa, the gyro unbalance term 
2 

= ACa , the gyro compliance term 
AP 


^9. “ 


(dimensionless) 


or (1/A#) At 
A(l/A#) 

(1/A#) 


the gyro scale factor term 


= A(G^» B^), the gyro alignment term 
= AQo:, the gyro nonlinearity term 

3) a = g ^ 

A#= 12 (sec/pulse) 

= 43200 (deg/hr) 

= 15 (deg/hr) 

5 2(l-cos 27rf At) ( „ j 

4) Pp(t> = j P„(f).cos2,„P^3(f). (C0^.ine„)2pj (Oj 

9 

= power spectral density of AP ((deg/hr) /cps) 

5) |Ts(3w)|2 = 1 

6) cos 0Q= 0, gyro horizontal 

= 1, gyro vertical 

sin 0 ^ = 1 /a [27 gyro horizontal or vertical 






CHART 2-4 


ACCELEROMETER COEFFICIENT ERRORS FROM STATISTICAL NOISE 


2) (^4^rms “ 


= (dimensionless) 

(gD^)At 


= , the accelerometer scale factor term 

= A{A.* Bj), the accelerometer alignment term 

= ADga, the accelerometer second order term 

= ADga^, the accelerometer third order term 
AP 

(g's) 

^ D^At 

= ADq, the accelerometer bias term 
a = g 

= 254 ((pulses/sec)/g) 

I 12 2(l-cos2ff f At) ( 9 

Pp(f)= ——3 9o> 

p 

= power spectral density of AP (g^cps) 

Tg(]cr)|2 = 1 

sin 9 q = 1 , for accelerometer horizontal 
= 0, for accelerometer vertical 






At (mm) 



Figure 2-1. Gyro Bias Error vs Time 
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^2 


rA(6; *B. ) = AQ<u 

(Va^) ' 



Figure 2-2. Gyro Scale Factor Error vs Time 
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E 

• Geometrical assumptions: g defines the vertical and earth rate w is a constant 
vector in inertial space, at an angle of 45° to the vertical. 

• The gyro is fixed relative to the turntable, which rotates at a constant angular 
velocity ^ relative to the laboratory. 


2. 2. 4 Bubble Level Compensation 


In the course of the computations required for Figures 2-3 and 2-4, it was noticed that, 
for the horizontal position of the accelerometer, the dominant term was the pickup of a 
variable component of g due to rotational environment noise. This effect could be re- 
duced if a bubble level were used periodically to measure the low frequency rotational 
motion during the calibration. Correcting for this motion, either mechamcally or 
mathematically, would reduce the low frequency part of the rotational noise spectrum. 
Two models were tried: 


1) The rotational noise spectrum, P 0 (f), was reduced to zero below a frequency 
corresponding to a 50- minute period. This result is tagged by the following 
symbol in graphs of accelerometer coefficient errors J 33 and in Figures 2-5 
and 2-6: I 



2) The rotational noise spectrum, P 0 (f), was assumed to be the squared modulus 

of a first order transfer function having an rms noise in A 8 of 4. 5 seconds of 

-2 

arc and a half-power frequency of 10 cps. This result is tagged by the 
following symbol in the figures: 



The results under these assumptions are plotted in Figures 2-5 and 2-6. 


2.2.5 Concluding Remarks 

Care should be exercised in drawing quantitative conclusions from the curves of noise- 
induced calibration errors. The noise data on which these curves are based are too 
scarce and too scattered geographically to be considered as the noise environment 
during a real calibration test. However, the curves can be used to support qualitative 
or comparative inferences such as the following: 

• A longer averaging time will reduce the calibration errors due to environment 
noise. 

• Calibration accuracy is strongly dependent on sensor orientation for both gyros 
and accelerometers. 
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• The large errors found for the horizontal accelerometer calibration can be re- 
duced considerably with bubble level corrections before and during the test. 

• To the extent that the assumed power spectra represent the actual test environ- 
ment, these results support Weinstock's conclusion* that the test bed should be 
isolated from the rotational (tilt) noise environment if a relative accuracy of 
the order of 10"6 is desired in the calibration of the sensor coefficients, 

• Subject to how realistic these power spectra are, it may be concluded that, for 
the gyro coefficients, and for the vertical accelerometer, the quantization 
error predominates. 


*H. Weinstock, loc, cit. , p. 45. 
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SECTION 3 



SECTIONS 

INTRODUCTION TO ALIGNMENT TRADE-OFF STUDIES 


This section serves as an introduction to the alignment trade-offs described in subsequent 
sections, hi the first subsection, the trade-off analyses are partitioned into three parts, 
each part corresponding to the contents of one of the Sections 4, 5, and 6. hi the second 
subsection the four nominal ISU orientations which are used in all alignment analyses are 
introduced. The last subsection describes the worst-case alignment errors due to instru- 
ment quantization. 

3. 1 ALIGNMENT TRADE-OFFS DEFINED 

The introduction to this document (Section 1) describes the alignment trade-offs which are 
accomplished in this study. The tabulation of the trade-off categories in the introduction 
is dictated by the listing as it exists in the contract statement of work, hi accomplishing 
these trade-offs, it is very convenient to partition the presentation into three parts. The 
partitioning is dictated by the functional descriptions of alignment as shown in Chart 3-1. 

The contents of Chart 3-1 are described in detail in Section 3 of the Development Document. 

For the purpose of this section, a brief description of the three separate routines is 

repeated. The three diagrams correspond to the three alignment techniques under study. 

Each diagram contains a routine entitled "Alignment Matrix Computations”, This routine 

corresponds to the mathematics which evaluates the elements of the alignment matrix as a 

function of the indicated inputs. The "Preprocessing Computations" transform the inertial 

instrument outputs into integrals of the measured values of applied acceleration or both 

applied acceleration and angular velocity components. The "Estimation Routine" produces 

estimates of £ and w from the outputs of the preprocessing routine. The coefficients 

of the ’Estimation Routine" are set by the "Estimation Matrix Computations" to provide 

E 

optimum estimation of £ and w from inputs corrupted by both environmental and 
quantization noise. 

The alignment trade-off problem is basically a question of determining alignment accuracy 
and time as functions of certain sources of error in the alignment technique. Thus there 
will in general be three cases corresponding to the three alignment techniques indicated 
in Chart 3-1. Moreover, since each technique can be considered as a collection of separate 
routines, a further breakdown of the alignment trade-off study by consideration of the 
several separate routines will be convenient. Each of the error sources to be considered is 
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either an error in the input to one of the routines or an approximation in the arithmetic 
employed in a routine. More explicitly. Mirror Alignment accuracy is a function of; 

1. Autocollimator readout errors. 

In Level Alignment both accuracy and time are functions of 1. plus; 

2. Environment (acceleration) noise 

3. Accelerometer readout quantization 

4. Accelerometer calibration accuracy 

5. Estimation Technique 

6. Estimation Computation (word length). 

In Gyrocompass Alignment both accuracy and time are functions of 2. through 6. above, 
plus. 

7. Environment (rotational) noise 

8. Gyro internal noise 

9. Gyro readout quantization 
lOo Gyro calibration accuracy. 

The trade-off analysis in this study is thus directed toward discermng the effect of these 
error sources on the accuracy of the alignment matrix obtained. However, autocollimator 
errors are not covered in this study because of unavailability of laboratory autocollimator 
data; and therefore Mirror Alignment will not be discussed in the succeeding sections. 

We see from Chart 3-1 that in the last two cases T is explicitly a function of g, and in 

E ^ 

the last case, o) . The trade-off study begins in Section 4 by relating errors in T to 

T? 

errors in estimating g and w . Sections 5 and 6 then relate the errors in the estimates 
E — — 

of g and w to the alignment error sources 2. through 10. listed above as they apply. 
Section 5 considers all of the error sources except inaccuracies in calibration constants, 
which are considered in Section 6, By combimng the results of Section 4 with the re- 
sults of Section 5 and 6 a complete picture of alignment accuracy and time, as a function 
of the error sources 2. through 10., is obtained. 

3. 2 ALIGNMENT ORIENTATION (CASES 1 TO 4) 

As will be seen in the general error equations of Charts 4-5 and 4-6 (Section 4), the 
appearance of terms such as (U* ^), (W* ^) etc. illustrates the dependence of align- 
ment accuracy on the nominal orientation of the body axes relative to the earth. 
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Experience has shown that four specific orientations bracket the extremes of this depen- 
dence. These four orientations are defined in Charts 3-2 through 3-5 which follow. All 
al:^nment accuracy analyses are developed for all four of these orientations. 

3. 3 AUGNMENT TIME VEESUS ALIGNMENT ACCURACY FOR WORST- CASE 
QUANTIZATION 

In Chart 3-6 we develop the relationship between worst-case quantization errors and 
errors in the estimate of g and as a function of sample time. Alignment errors in 
terms of errors in the T matrix, AT, are obtained by substitution of these worst- 
case quantization error expressions into the error matrices of Section 4. 2. Alignment 
errors in terms of total rotation ai^le or cone angles are obtained by substitution of the 
elements of the ATT*^ of Section 4. 2 into the expressions developed in Chart 4-3. 

Application of these results show that the gyrocompass azimuth error, due to a worst- 
case, one pulse, "east gyro" quantization error (A^ = 12 seconds of arc per pulse) is 
on the order of 3700/ At seconds of arc ( At is in mmutes). On the other hand, the level 
error due to worse-case, one pulse, accelerometer quantization error (D^ = 254 pulses 
per g) is on the order of 37/At seconds of arch (At in minutes). Thus the gyrocompass 
alignment error is dominated by "east gyro" quantization. An alignment time on the 
order of two hours is therefore required to reduce the azimuth error to the order of the 
level error. 

Because of this gyro quantization effect, the alignment studies of Section 5 are devoted 
principally to an investigation of noise filtering techniques for level alignment. 
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CHART 3-2 


NOMINAL ALIGNMENT ORIENTATION 



cJE . 

Bi = 0 

£* §1 


Bg = o)Esm X 


cijE • 

cos X 

g’ ?3 


when X - 45° 

= 0 

Bg == 0.707 
££^•63 = 0.707 

£■§1 = 0 

I * §2 = 0 

I* §3 = e 
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CHART 3-3 



3-6 




CHART 3-4 



3 



CHART 3-5 


NOMINAL ALIGNMENT ORIENTATION 



kflTs cosX 


[T] = 


\l2/3 
kUTi sinX 




E 


Bj =>H73 

Bg 

Bg = -vjlTs 


,E 


when X = 45° 


Case 4 


Bg and Bg are equal angles from earth rate 
Bj^ IS in the earth rate- east plane 


\l 1/3 cosX - \ll/2 sinX \Jl/3 cosX - \jl/2 sinX 

-xfiTe -nUA 

xTiTs sinX + '^1/2 cosX \[i73 sinX - \il/2 COSX 


g • = (-vjl/B COS X) g 

g • Bg = (>|l/3 cos X - >11/2 sin X) g 
g • Bg = (>jl/3 cosX +>1 1/2 sinX) g 



4176 

(4176 - 1/2) 

(41/6 + 1/2) 


0. 407 

-0. 093 

0.907 

[T] = 

vf^ 

-n|i76 

-nUTs 

= 

0, 815 

-0. 407 

-0. 407 


4i/6 

(Nfi/e + 1/2) 

(4 1/6 - 1/2) 


0. 417 

0.907 

-0.093 


tc®- B^ = 0.577 w® 
• Bg = 0. 577 a® 
Bg = 0,577 


g- Bj = 0.407 g 

g • Bg = -0. 093 g 
g • Bg = 0. 907 g 
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CHART 3-6 


QUANTIZATION ERROR 

Assumptions 

• Linear instruments 

• No noise 

• Perfect calibration 

•V. 

• Constant gravity and earth rate inputs 

• Instrument and body axes are perfectly aligned 

These assumptions say that 

£ . B = (l/D^)(P^/At) 

Assuming a n^ quantum error in and a n^ quantum error in P^ we have 

coE o^^At 


A(g • Bj^) 


g D^gAt 

Substituting the nominal scale factors and 


C' A 

n = 1 pulse n = 2 pulses 


we have 


A{^ . Bj^) 

133 X 10"'^ 

o>E 

At(in minutes) 

A(g • Bj^) 

13 lx 10"® 

g 

At(in minutes) 
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SECTJON 4 



SECTION 4 

ALIGNMENT ACCURACY VERSUS g, PRECISION 


We begin this section by defining a basic measure of alignment errors (the elements of 
the Att^ matrix). We then relate this basic measure, through the T matrix, to errors 
in estimates g and w® for each of the four orientation cases and for both level and 
gyrocompass alignment. Finally, we develop the statistics of the basic measure of the 
alignment errors in terms of the statistics of the errors in the estimates of g and 

The results of this section will be used in Section 5 and 6 to transform errors in g and 
E 

« to equivalent errors in the defined basic measure. 

4. 1 GENERALIZED ERROR EQUATIONS 

In the Development Document, alignment was defined as the initialization of the T matrix 
which transforms from an ISU fixed set of axes to an earth-fixed, local level set. Three 
types of alignment techmques were presented: Mirror Alignment, Level Alignment and 
Gyrocompass Alignment. Errors m Mirror Alignment are directly a function of the 
autocollimator survey and not a function of the outputs of the inertial instruments or 
computations. For this reason, the following generalized error equations are developed 
for Level and Gyrocompass Alignment only. 

The T matrices for level alignment and gyrocompass alignment are reproduced for refer- 
ence in Charts 4-1 and 4-2, respectively. It should be noted that the computed matrix 
IS always orthonormal irrespective of the errors m U and W. Thus, the alignment error 
can be e 3 q>ressed as the difference between two sets of orthogonal axes — the true and 
computed sets. This difference can be esq)ressed by a single rotation which aligns the 
two sets. This rotation will be called the "total rotation angle" in the discussions which 
follow. The difference can also be expressed by the three angles between pairs of 
corresponding axes in the two sets. These angles will be called "cone angles". 

Referring to Chart 4-3, the cone angles are approximated by the magnitude of the cross 
product between corresponding axes in the true and computed sets. The equations for 
the squares of the cone angles given in Chart 4-3 are developed in Chart 4-4. 

T 

Both the cone angles and the total rotational angle are developed as functions of ^TT , 

As the errors in g and w will be small, the AT matrix is obtained by taking the first 
difference of the T matrix in terms of g and to . The At matrices for level alignment 
and gyrocompass alignment are presented in Chart 4-5 and 4-6, respectively. 
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CHART 4-1 


LEVEL ALIGNMENT MATRIX 
Inputs (g. B^'), (g“ Bg), (g* Bg) and 


From these quantities the alignment matrix is given by: 



where 

• (Mi*U) = (U-B^) 

« iI^xUI= Cl - 

• (y. = (g- Bj^)/g 

e g = C (g- + (g- Bg)^ + (g . Bg)2 ]l/2 

An optional technique might utilize any of the following additional inputs: 

• The zenith angle (6^) of mirror one might be utilized to find (Mj^ • U) from 

(l^* U) = cos 0^ 

• The magmtude of gravity (g) might be supplied from a local survey. This piece 
of information can be utilized to reduce the number of required accelerometers 
to two. 
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CHART 4-2 


GYROCOMPASS MATRIX 


Inputs (go (g' ^ 2 )^ (g* Bg)^ (w®. Bj^) (w^< Bg)^ and Bg) 
From these quantities the alignment matrix is given by: 




.E 




0 

1 

0 


(W-Bj) 

(W-Bg) 

(W-Bg) 

T 


0 

0 

1 

Iwxul 


(U-Bi) 

(U-Bg) 

(U-Bg) 



1 

IWxUl 

(W'U) 

Iwxul 

0 


(W X U)*(B2X 

Bg) (WxU)-(B^xB 

(WxU).(BjXB2) 


where 

(W . U) = (W . B^)(U. B^) + (W • B2)(U • Bg) + (W • Bg)(U. Bg) 

|WxU| = Cl - (W-U)2]1/2 
(W.^) ={w®.Bj^)/oF 
(U.^) =(g.Bj^)/g 

tsP - C(w®" B^)^ + (w^- Bg)^ + (w®. Bg)2]^/2 
g -C(g.B/+(g-B2)2+(g.B3)2]l/2 

An optional technique might utilize any of the following additional inputs: 

« The local latitude (X) might be utilized to find (W * U) from 

(W • U) = cos X 

e The magmtude of gravity (g) might be supplied from a local survey. 

• The magnitude of earth rate (w^) might be supplied from a local survey. 


A use of all additional inputs could reduce the number of necessary instruments to three 
(either two accelerometers and one gyro, or one accelerometer and two gyros). 
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CHART 4-3 


ALIGNMEN T PRECISION 

The alignment matrix has been defined as a transformation from body to earth axes- 
That is 


El 


U 





^2 


E 

= 

T 



1 

IM 

1 


N 




?3 


If the elements of T are in error, an erroneous earth frame will be defined, or 






El 


= 

T + AT 



l3 




?3 


A multiplication of this matrix by T^ will yield a transformation from the real earth 
axes to the erroneous axes. 


pi1 










ll 

E' 

-2 

= 

T -(-AT 


tT 


-2 

= 

I -f- att"^ 


?2 

E' 

-3 






_-3_ 






From Chart 4-4, which follows, the cone angle errors of up, east, and north are given by 


U X U'j ^ = 

att^ 

2 

12 + 

att"^ 

2 

13 

[ . 1 2 

r Ti 


r t1 

2 

[E’xEp = 

ATTpj 

I 2 I ^ 

ATTJ 

23 

In’xnP = 


2 ^ 1 
31 + i 

Iatt'^J 

12 

I 32 


In all error analyses in this document we have constrained the computed (E^) frame to 
always be orthonormal. The real (E-) frame is by definition orthonormal. Therefore, 

rp ^ 

the ATT matrix (to first order) will always be an anti^mmetric rotation matrix. The 
three independent quantities in that matrix will therefore represent the up, east and north 
components of the small angle rotation vector. A representation of the rigid bocfy rotation 
between the computed and real frame will therefore be 

(total rotation angle)^ = (ATT^)j 2 + (ATT'^)jg (ATt ‘^)23 
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CHART 4-4 


CONE ANGLES FROM A TT^ MATRIX 


U X U’l ^ = (U X U') • (U X U’) 

= (U* • U’) - (U • U')^ 


From Chart 4-3, 


= I + ATT-^ E 


From which it follows that 


y' - U = (ATt'^)jjU + (ATt”^)j2E + (ATt'^)^ 3N 
(U' - U) • (U' - U) = (U’ * U’) - 2(U' • U) -r 1 


= (ATt'^)jI ^ + (ATt'^)j2 ^ + (ATt'^)i3 


I^U’ - U) ■ U 2 = (U' . U)2 - 2(U' • U) + 1 


From which it follows after substitution 


UxU’ ^ = (ATt’^)j^^ + |TaTt'^)j 3 ' 


In the same way, 


T) 2 
^23 


lExE'^ = |^TT^)2j ^ + ( ATt '^)23 
NxN'l^ = |^TT^)3^ ^ + (ATT*^)32 
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4-6 


LEVEL ALIGNMENT ERROR EQUATIONS 


Assuming • (g* B|), (g* Bg), and (g* Bg) are m error 

• No constraints, such as the magnitude of gravity, are used in the solution of T 

Then, the first order error in the alignment matrix is given by 


LAT] = 


10 0 


0 sin a, cos 


ctj^ cos 


0 -cosoj^ sinoj 


A(U.Bi) 




A(U-Bj) - 


(U-B^){U-B3) A(U-B3)' 

A(U-Bj^)+ 

L lUxB^lJ 

(U.Bg) (U-B^) 

A(U* B^) + A(U- Bg) 


A(U. I 


(U-Bi)(U.Bg) 


-A(U- 


L 




UxB, 


'(U-Bg) 


UxB^ 


A(U.Bj) 


where A(U-B(^) = S L 5^^ - (U* B^)(U • 1 


■A(g.Bj 


, k = 1, 2, 3 


and 


(U* for (k = l, 2, 3) are the zenith angles of the body axes. 
The rotation matrix IS given by [aTt'^] = [at] It*^] 


A{U Bg) 

B j^) + 

i!!=‘5iL 

(U- B,) 
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CHART 4-6 




4. 2 ERROR EQUATIONS FOR CASES 1 TO 4 


In this subsection the error equations are evaluated for each of the four selected orien- 
tations. In the eight charts (4-7 throi^h 4-14) which follow, the At and Att"^ matrices 
are listed for the four orientations for both level and gyrocompass aligimient. 


E 


'4.3 STATISTICS OF T FROM THE STATISTICS OF g AND w 

The one -sigma of the elements of T and Att^ are derived from the statistics of 
E 

g and ^ as follows. 

rpi 

An element of either AT or Att can be e^qiressed in the form 

A Ck^(£-Bk> • gk ) } 

I B " ' 


k 


where Cj^ and Cj^^g are constants. The one-sigma value of this ejqpression is 




a = \ E 




^1 






- + 


- fi 


CO 


E 


k+3 


where E [a] is the expected value of A, 

’ Ik) 


^k “ ® 


^k+3 


= E 


• 5k)' 


0 


E 


As before, we assume that the functions 

A(g • Bj^) 


A(^^ • Bj^) 


and 


g 


U) 


£ 
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LEVEL ALIGNMENT 




CASE 1 



A(g.B^) 

^(£• 12 ) 




1. 000 

1. 000 

0.000 



g 

g 






A(g.Bj) 


AT = 

0. 000 

0. 000 

-1. 000 





g 





^(£•§ 2 ) 



0. 000 

0.000 

-1.000 



— 


g 



- 

A(g.B^) 




0 non 

1 non 

1 nno 




g 

g 



A(g.Bj) 




att-^ = 

-1. 000 

0.000 

0, 000 



g 





ACg.Bg) 





-1. 000 

0,000 

0. 000 



g 
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CHART 4-10 


LEVEL ALIGNMENT 


A(g.B2) 

0,833 +0.037 0.371 

g e g 


AT = 


A(£-B3^) A(£.B2) A(g.Bj) 

.0.333 0.015 +0, 148 

g g g 


A(£.B^) A(£.B2> A(£.Bg) 

.0. 167 -0. 007 +0. 074 


CASE 4 


A(£'Bi) A(£.Ba) A(g.^) 


A(£-Bi) A(g.B2) A(g.^) 


0.037 +0.992 

g g 


+0,083 


A(£.B^) A(£.B2) A(£.B3> 

0.075 0.442 0.078 

g g g 


-0.371 +0.083 

g 


- +0. 175 


g 


A(£.B2> A(g.B3> 

0. 037 0. 098 0, 027 

g g g 


A(£.Bj) Atg.Bg) A(£-Bg) 

-0.742 +0.412 +0.375 

g g g 


A(£.B^) A(£.B2) A(£.^) 

-0.371 1.007 +0.065 

' g g g 


0.000 


ATT 


,T _ 


A(s-Bj) A(£.B2) A(£-B3) 

-0.816 +0.408 +0.408 

g g g 


A(£.Bj) ACg.Bg) A(£.B3) 

-0.408 0.908 +0.092 

g g g 


A(£-Bi) ACg.Bj) A(£.B3> 

0. 816 0. 408 0. 408 

g g g 


0. 000 

A(£.B2) A(g,^) 

0. 445 +0, 045 


A(g-Bj) A(£.B2> A(£.B3) 

0.408 +0.908 0.092 

g g g 

A(g.B2) A(g.B3) 

-0, 445 0. 045 

g g 


0.000 
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GYRO COMPASS 



CHART 4-11 
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CHART 4-13 











CHAIIT 4-! 4 


AT = 


ATt"^ a 


GYRO COMPASS 


CASE 4 


0 833 +0 037 -0 371 

e s 


A(w^Bi) Afw^Bg) ACu^Bg) 
-0 471 — +0 236 +0 23$ ~ 


Bl> At£.^> 

-0 833 0.037 +0 371 

I! e g 

A(u^Bj) 

^ 943 0.471 0 471 — H— 

W® 


A(£ Bj) A(j^ Bg) A(g 

0 037 — — *0 992 *0 083 

g g g 


A(ijDi) AljrBgJ A(t^) 

0 816 0 408 0 408 

g g t. 

A(w^Bi) Alu^Bg) 

-1 049 ^ tO 524 =- +0.524 =— 

A(£-Bi) A(j£ Bg) 

0 371 0 083 0 175 

g g 8 

A(yEBp A(u^B 2) A{y^Bg) 

-0 471 ~ +0 236 +0 236 

,,E E E 


A(^ Bj) A<j, 

-0.371 rO 083 -0 175 

g g g 


A(g.B,) Aig Bj) A(^ Bg) 
-0.816 -»0 408 0 408 

g g g 

40 106 =— -0 053 =— -0 033 — - ■ — 

S E E 

M w w 

_ Atg.Bg) 

-0.037 -0.992 — -0 083 

g g g 

A(y^B|) 

-0 471 tO 236 -0 236 

SEE 

(I) u u 


A(gB^) A(g^) A(g^) 

0 816 0 408 0 408 


A(b B j) A(£ A(^ 

0.408 ■ ■ . -0 908 -0 092 — 

U g g 


A(g B.) A(£.^) 

-0 816 1-0,408 0 408 

g g g 


i A(£ B,) Al£ A(£ Bg) 

! 0 816 0 408 0 408 


AU^B,) Afuj^Bg) 


A(£-Bi) A(£.^) A<£ ^ * 

, -0.816 0.408 -0 408 

A(g.y A(£ Bg) A(£.Dg) g g g 

-0 408 -0 908 +0.092 I c 

g g g I a(u,®Bi) §3) 

i *l I55 —-0 577 0.577 ! 

r ^.2 B B * 

I <A, Cl. I 


4-16 






have a zero mean and are uncorrelated so that, 
finally: = 0 
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SECTION 5 



SECTIONS 

ALIGNMENT PROCESSING TECHNIQUES 


In this section we investigate the effects of noise, quantization and computer word length 
upon alignment accuracy. The analysis is based on a Monte Carlo simulation which is 
described in Section 5,1. In Section 5.2 the general properties of several estimation 
techniques are developed and a "recommended technique" is selected. These techniques 
are investigated in detail in Sections 5. 3 and 5. 4. In particular, we investigate the 
following relationships: 

1) Alignment accuracy and alignment time versus sample rate 

2) Alignment accuracy versus number of iterative steps 

3) Alignment accuracy versus sensor quantization 

4) Alignment accuracy versus computer word length. 

We briefly summarize the results of these investigations in the following paragraphs. 

The recommended technique for level al^nment is the posterior -mean estimate of the 
instantaneous components. An alternate technique is the simple average. These 
recommendations are based on the results presented in Tables 5-1 and 5-2. The 
posterior -mean estimate is not sensitive to the distribution of the environment noise, 
i. e. , gaussian or nongaussian. 

Rotational motion from the environment is the dominant source of level alignment error 
for long integration intervals of the order of one minute or greater. Quantization and 
sensor noise are dominant for short integration intervals, less than one -fourth minute. 
The posterior -mean estimate of the instantaneous gravity components is more accurate 
than the simple average m most cases of interest. For very short intervals (on the 
order of 15 seconds) these methods have comparable accuracy. The posterior -mean 
estimate is less accurate than the simple average for large quantization and small 
mtegration intervals, the order to 30 seconds. 

For Gyrocompass, the simulation results confirm the conclusions of Section 3.3 — that 
the alignment error is dominated by gyro readout quantization. For this reason, only 
qualitative relationships between al^nment accuracy, alignment time, sample rate, 
sensor quantization and computer word length are developed. 
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These results employ the environment data given by Welnstock. * It is recommended that 
the motions of a proposed test laboratory be studied, since deviations from the nominal 
environment described by Weinstock will change the characteristics of the recommended 
techniques. Deviations from the nominal environment will be most likely in magnitude 
of the spectrum and not in the shape. In other words, low frequency motion is larger than 
the high frequency motion. The magnitude cf the low frequency motion depends on the 
location of the test bed, the time of day, local human activities, etc, it may change by 
an order of magnitude. 

5. 1 DESCRIPTION OF SIMULATION PROGRAM 

The function of the computer simulation is to generate data similar to that which is 
obtained from the accelerometers and gyros, to process the data with various estimation 
techniques, and to compare the true alignment matrix with the estimated alignment 
matrix. A functional description of the simulation appears in F%ure 5-1. A detailed 
description of the simulation appears in the report describing the LABSIM Program. ** 

A functional description of the simulation is presented in Figure 5-1. The components 

F 

of gravity, g , and earth rate, w , in the earth frame are transformed to components 
in the level frame, i. e. , a frame nominally aliped with the earth frame but moving 
with the laboratory. This transformation depends on the low frequency (LF) rotational 
motion of the laboratory (10" to 10" cps). The high frequency (HF) motion ( 10“ to 
10 cps) does not significantly change the orientation of the level frame, because the 
amplitude is smaller (see Figure 5-2). The resultant components, g* L-(t) and 
« • L- (t) , are integrated to simulate the integral readout of the sensors. (In the 
simul^ion program, the integration is done analytically. The integrals serve as the 
input to the program.) The integral of the HF rotational noise is added to these 
components, and the sum is transformed to the body frame. In the simulation we 
assume that the sensor input axes are parallel to the body axes. Gyro noise is added 
to the gyro components. No accelerometer noise is added. The resultant signal and 
noise is quantized. (Quantization is accomplished by dividing by the appropriate scale 
factor, rounding and then multiplying by the same factor.) The outputs of the quanti- 
zation routine simulates the actual sensor ouiputs. 

The estimation routine computes the estimates of gravity and earth rate in the body 
frame, g* and w • B^. The estimation routine uses certain estimation matrices 


*H. Weinstock, "Limitations of Inertial Sensor Testing Produced by Test Platform 
Vibrations", NASA Electronics Research Center, Cambridge, NASA TN D-3683, 1966. 

**The LABSIM program is described in "Simulation Program of Inertial Sensing Unit 
for Evaluation of Alignment Techniques", a technical report prepared by Univac, 
Aerospace Analysis Department, St. Paul, Minnesota, January 1968. 
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Figure 5-lo Functional Description of Simulation 











M and b. The matrices are computed as part of the simulation, but in a laboratory 
system M and b will be computed before the alignment. The quality cf the estimate 
IS measured by the alignment total rotational error. This error is defmed m Section 4 
as the magnitude of the su^le rotation required to bring the estimated earth frame into 
coincidence with the true earth frame. 

In the following subsections, various routines m the simulations are discussed: test 
inputs (5. 1. 1), estimation routme (5.1.2), estimation-matrix routine(5. 1. 3), rotational 
alignment error (5. 1.4). 

5. 1. 1 Test Inputs 

The laboratory environment introduces tivo types of "noise" — rotational and translational 
motion. The acceleration introduced by the translational motion was not simulated, 
since it is much smaller than the accelration modulation due to rotational motion. 

The rotational motion was divided into two components — high frequency and low frequency. 
The high frequency component was formed with a random-number generator which had 
a gaussian distribution. The power spectrum was shaped with a recursive filter to give 
the spectrum illustrated in Figure 5-2. Independent HF motions were applied to the 
N-S axis and the E-W axis. The rms amplitudes were 0, 8 second cf arc. A nongaussian 
high frequency rotational noise was also simulated. A bimodal density function was used 
for this purpose with the second peak containing one tenth of the total probability. 

The nongaussian random -number sequence was formed from the original gaussian sequence 
by replacing every tenth number on the average with a gaussian variate whose mean is 
1. 6 seconds of arc and variance is 0.8 second of arc. The variance of the nongaussian 
noise was greater than the variance of the gaussian noise. 

The low frequency rotational components were simulated with a harmonic motion. Two 
motions were used: 1) one hour period and one mmute of arc amplitude, and 2) one-half 
hour period and 30 seconds of arc amplitude (see Figure 5-2). The axis of rotation was 
chosen in the horizontal plane, 45° north and east. The experimental data in Figure 5-2 
was obtained from Weinstock.* 

The rms noise level of the gyro is 0. 005°/hr, It was formed with a random -number 
generator which had a gaussian distribution. The power spectrum was shaped with a 
recursive filter to give a half -power frequency of 2 x 10“® cps and a roll-off of 
6 db/Oct. The noise applied to each axis was independent of that on the other two axes. 


*H. Weinstock, "Limitations on Inertial Sensor Testing Produced by Test Platform 
Vibrations", NASA Electronics Research Center, Cambridge, NASA TN D-3683, 1966. 
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ANGULAR DISPLACEMENT POWER SPECTRUM 


ANGULAR VIBRATION ENVIRONMENT 



FREQUENCY (CPS) 


Figure 5-2. Simulation Rotational Inputs 
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5.1.2 Estimation Routine 


The estimation routine implements the basic eq,uation for the estimated components, e.g. , 
in level alignment 


= MX + b (5-1) 


where X is the preprocessed measurements 
At KAt 

J (a.B.)dt, • • • , J (a*B.)dt; i = 1,2,3 
■ 0 (K-l)At 

The matrices M and b are evaluated in the M and b Evaluation Routine. The computer 
word length restriction is mtroduced in the evaluation of (5-1). Word lengths of 27, 24 
and 15 bits were used. A detailed discussion of the estimation routme is presented m 
the report on LABSIM. 

5.1.3 M and b Evaluation Routine 


The elements of M and b of equation (5-1) are computed m the M and b Evaluation 
Routine. These matrices are m most cases functions of the number of samples K, the 
inter sample time At, prior measurement of the noise power spectra, initial estimate 
of the alignment matrix, and prior measurements of gravity and earth rate. The 
equations for M and b depend on the alignment technique being used. The simple 
average estimation techniques do not require all of the above quantities to evaluate 
M and b. The basic estimation equations are developed in Section 5 of the Development 
Document. 

5. 1. 4 Error Evaluation Routine 


The Error Evaluation Routme computes the rotational alignment error defined in 
Chart 4-3. This angular error is the rotation required to bring the estimated earth 
frame into coincidence with the true earth frame. The absolute magnitude of this ro- 
tational error will be used to compare various techniques. Note that the earth frame 
is moving relative to the body frame. The alignment error is based on the predicted 
orientation of the earth frame one second after the last measurement. This is a 


S* 

/\ 
I- §2 
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worst-case delay for initializing the navigation program. The errors presented in the 
following subsections are the rms rotational errors based on 10 independent Monte Carlo 
trials. 

5.2 SELECTION OF RECOMMENDED TECHNIQUES 

In Section 5 of the Development Document, several alignment estimation techniques were 
derived. We will select a "recommended" estimation technique in this subsection, for 
both Level and Gyrocompass Alignment, from among those derived in the Development 
Document. The selection is based upon three criteria; 

• The sensitivity of alignment accuracy to the noise distribution, gaussian and 
nongaussian. Only techniques which are not sensitive to the noise distribution 
will be considered further. 

• The relative accuracy of the techniques for different integration times, ^t; 
orientations, TEBj and number of samples, K. 

• The general computational requirements which include complexity, accuracy 
(double precision, etc.) and setiq) procedures for laboratory test. 

The selected techniques can be considered "recommended" only to the extent that the 
actual laboratory environment noise approximates the simulated noise. (The noise 
simulation was based upon the environment data given by Weinstock, loc. cit. ) Deviation 
from this nominal environment could change the characteristics of the estimation 
techniques enough to change our choice for recommended technique. Specifically, if 
the low-frequency motion were small (on the order to five seconds of arc), then our 
choice would be different. 

The following parameters were used in the test simulations. Nominal quantization was 
introduced: 1.27 x 10" ft/sec and 1. 22 x 10" rad. No word-length restriction was 
used in the estimation routine. The low-frequency rotational motion was one minute of 
arc with a one hour period. Inputs to the Estimation Routine were taken symmetrically 
about the zero phase of the LF motion, i.e. , maximum angular velocity. A nongaussian, 
high-frequency, rotational noise was also introduced. It had a bimodal density function 
with the second peak containing one tenth of the "total probability" (see subsection 5. 1. 1). 
Two orientation cases are considered — 1 and 2 (see orientation cases in Section 3). 

5.2.1 Level Alignment 

In the Development Document four estimation techniques for Level Alignment were 
presented: simple average, least-squares and posterior -mean estimate of average 
components, and posterior -mean estimate of instantaneous components. The least-squares 
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estimate has not been investigated in detail, because it is a special case of the posterior- 
mean estimate, from a functional viewpoint (see Section 5 of the Development Document). 

The simulation results are presented in Tables 5-1 and 5-2 for gaussian and nongaussian 
HF noise, respectively. The table entries are the rms rotational alignment errors ex- 
pressed in seconds of arc for 10 trials. In all cases the instantaneous estimates are 
superior to the estimates of the average components. The accuracy of the instantaneous 
estimate is not sensitive to the HF noise distribution. 

The posterior -mean estimate of the instantaneous components is selected as the recom- 
mended estimation technique for level alignment. The simple average is an alternate 
technique because of its computational simplicity. The posterior -mean estimate can be 
used iteratively. The characteristics of this technique are presented in Section 5. 3. 2. 

5.2.2 Gyrocompass Alignment 

In the Development Document three estimation techniques for Gyrocompass were 
presented: simple average, least-squares estimate of average components, and 
posterior -mean estimate of average components. The least-squares estimate is not 
investigated in detail, because it is a special case of the posterior -mean estimate, from 
a functional viewpoint. 

The simulation results are presented m Table 5-1 and 5-2 for gaussian and nongaussian 
noise, respectively. The table entries are the rms rotational alignment errors expressed 
in seconds of arc. The posterior -mean estimate is not significantly better than the 
simple -average estimate. Since the dommant source of error in Gyrocompass Align- 
ment is gyro quantization (see Section 3. 3), this result is expected. The simple average 
is therefore selected as the recommended technique. Note that if the gyro quantization 
effect could be made small, the recommended technique could be improved upon by 
using the instantaneous estimate of gravity (given in Section 5.2.1) and the average 
estimate of earth rate. 

5. 3 CHARACTERISTICS OF RECOMMENDED TECHNIQUES 

In the precediiig subsection, the posterior -mean estimate of instantaneous components 
was selected for level alignment. In this subsection, we investigate the characteristics 
of this technique and the simple -average technique. 

For gyrocompass alignment, the simple -average estimation technique was selected; 
its characteristics are also investigated. Specifically, the relationships are developed 
between: 
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TABLE 5-1 

SELECTION OF RECOMMENDED TECHNIQUES 


Rotational Alignment Error (sec) 

Gaussian Noise Distribution 
A - Simple Average 

B - Posterior-Mean Estimate of Average Components 
C - Posterior- Me an Estimate of Instantaneous Components 
D - Simple Average 

E - Posterior- Mean Estimate of Average Components 


K = 1 
At = 10 min 


ORIENTATION 


LEVEL 

GYRO 

B 

C 

D 

E 

31.1 

31.1 

260 

320 

47.8 

34. 6 

297 

324 


K = 5 
At = 2 min 


ORIENTATION 


LEVEL 


B 


GYRO 


K = 5 
At = 5 min 


ORIENTATION 


LEVEL 


B 


I 

60.0 

60a 

n 

66.8 

45.4 


GYRO 
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TABLE 5-2 

SELECTION OF RECOMMENDED TECHNIQUES 


Rotational Alignment Error (sec) 

Non- Gaussian Noise Distribution 
A - Simple Average 

B - Posterior- Mean Estimate of Average Components 
C - Posterior- Mean Estimate of Instantaneous Components 
D - Simple Average 

E - Posterior- Mean Estimate of Average Components 


K = 1 
At = 10 min 


ORIENTATION 


LEVEL 


GYRO 


A 

B 

C 

D 

E 

30.8 

30.8 

30.8 

260 

320 

34.4 

48.9 

34.4 

410 

332 


K = 5 
At = 5 min 


ORIENTATION 


LEVEL 

GYRO 

A 

B 

C 

D 

E 

60 

60 

34 

60 

364 

68 

45 

40 

207 

495 
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• Alignment accuracy, alignment tune, and sample rate 
0 Alignment accuracy and sensor quantization 

* Alignment accuracy and computer word length 

« Alignment accuracy and number of iterative steps (level aligmnent only). 

These are the topics of the following subsections. Since gyro quantization is the dominate 
source of error, we obtain only qualitative relationships for gyrocompass alignment. 

5.3.1 Alignment Accuracy 

Alignment accuracy is measured by the rms rotational alignment error. The alignment 
time is the total measurement interval KAt. Note that the computation time is small rela- 
tive to KAt. The object is to relate alignment error, alignment time, and the intersampie 
time At. 

Several tests were performed with the following parameters. Nominal quantization was 
introduced, 1.27 x 10"^ ft/sec and 1.22 x 10"^ rad. No word length restriction was used in 
the estimation routine. The low-frequency motion was one minute of arc with a one-hour 
period. The mputs to the Estimation Routine were taken symmetrically about the zero 
phase of the LF motion, i.e. , maximum angular velocity. A gaussian high-frequency noise 
was introduced. Orientation cases 1 and 2 were considered. In the following subsections 
a noniterative and an iterative level alignment techniques are discussed. Subsequently, 
gyrocompass alignment accuracy is described. 

5.3. 1.1 Noniterative Level Alignment 

The alignment error for the posterior -mean estimate of the instantaneous components is 
given m Figure 5-3 as afunction of alignment time KAt for K = 1, 2, 5, 10. For K = 1, 
the error is the same as that for the simple average. The error peaks at 30 minutes 
smce the low-frequency rotational noise has a one-hour period and data is taken symmetri- 
cally about vertical. The rotational alignment error is based on the instantaneous orienta- 
tion at the end of the measurement interval. For K = 1 and At equal to one hour, the 
rotational motion would be eliminated completely. Note that there is a significant error 
reduction in using K = 5 instead of K = 1. The reduction in going from K = 5 to K = 10 is 
not significant. In the following discussions we will use K = 5. Since data was taken 
symmetrically about zero phase, the error appears to approach zero as KAt approaches 
zero. The sampling and quantization are "in phase". A detailed investigation is required 
for small KAt. 
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Posterior -Mean Estimate of Instantaneous Components 

Gaussian Noise Distribution 

Quantization: 1.27 x lO"^ ft. /sec/ and 
1.22 X 10-'^ rad. 

Orientation* Case 1 



ALIGNMENT TIME, KAt (MIN.) 


Figure 5-3. Level Alignment Error vs Alignment Time 
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There are two basic types of errors: prediction error and noise error. The first type of 
error occurs since we attempt to predict the position of the gravity vector at the end of 
the measurement interval. When the measurement interval KAt is decreased, the predic- 
tion error is decreased. The noise error results from the HF rotational motion and 
accelerometer quantization. When the measurement interval KAt is decreased, the noise 
error is increased. The optimum interval is one for which these two errors are equal. 

In the following subsection, we find that the optimum At is of the order of thirty seconds 
for K = 5 and for the noise spectra given at the beginning of this subsection. Since the 
choice of At is strongly dependent upon the noise spectra, the final selection must be based 
upon measurements of the actual noise spectra in the laboratory. 

The alignment error for the simple average is presented in Figure 5-4. For simple average, 
only one sample is taken, i.e. , K = 1. Note that alignment accuracy is not strongly depend- 
ent on orientation. 

5.3. 1.2 Iterative Level Alignment 

The posterior -mean estimate of the instantaneous gravity components can be used in an 
iterative alignment technique. Initially we liave a prior estimate of the gravity components 
and the associated covariance matrix (see Section 5.4. 1. 3 of the Development Document). 
These are inputs to the first step of the iteration. Using the posterior -mean estimate, 
we obtain refined values for the gravity components and covariance matrix. These refined 
values are inputs to the second iteration step. The procedure continues in this way. 

This procedure was tested with different intersample tunes At and different noise inputs. 
Specifically, K = 5 and At = 0.25, 0. 5, 1. 0, 2,0 minutes. The quantization levels used 
were 1.27 x lO'^ ft/sec and 1.22 x 10'^ rad. The word length was not restricted in the 
estimation routine. The high-frequency rotational noise was gaussian. The low-frequency 
noise had a period of 58 minutes and amplitude of one arc minute, and in a second case, a 
period of 29 minutes and an amplitude of 0. 5 arc minutes. The initial input to the Estima- 
tion Routine was taken when the body frame was level. Orientation Case I was used for 
all tests. 

The simulation results are presented in Figure 5-5, 5-6, 5-7, and 5-8 for the 58-minute 
period. The simple average is also plotted for comparison. Note that the simple average 
Is not iterative. The ordinate is the rotational alignment error in seconds of arc. * 


*These errors are the result of a smgle trial and not the rms error of ten trials as used 
above. 
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Simple Average 

Gaussian Noise Distribution 

Quantization; 1.27 x 10“'^ ft /sec. and 
1.22 X 10-4 rad. 


Orientations: Cases 1 and 2 



Figure 5-4. Level Alignment Error vs Alignment 
Time For Simple Average 
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ROTATIONAL ALIGNMENT ERROR - SEC 



Figure 5-5. Alignment Error — Iterative Technique 
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Environment Motion- 



TIME - MIN, 


Figure 5-6. Alignment Error - Iterative Technique 
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Environment Motion’ 
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Quantization: 
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Gaussian Noise 
Orientation: Case I 
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Figure 5-7. Alignment Error - Iterative Technique 
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At = 1/4 Min. K = 5 
Environment Motion- 
Period = 58 Min, 
Amplitude = 1 Mm. 
Quantization: 

1.27 X 10-2 Ft/Sec 
1.22 X 10"'^ Rad 
Gaussian Noise 
Orientation; Case I 
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Figure 5-8, Alignment Error - Iterative Technique 
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At = 2 Min. K = 5 
Environment Motion* 
Period = 29 Min. 
Amplitude = 30 §ec. 
Quantization* 

1.27 X 10-2 Ft/Sec 
1,22 X 10-4 Rad. 
Gaussian Noise 
Orientation: Case I 
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Figure 5-9. Alignment Error - Iterative Technique 
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Figure 5-10, Alignment Error - Iterative Technique 
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Figure 5-11. Alignment Error - Iterative Technique 




Figure 5-12. Alignment Error - Iterative Technique 
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Figure 5-13. Gyrocompass Alignment Error vs Alignment 
Time For Simple Average 
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The absissa is time in minutes. The curves represent the alignment error based on data 
taken over the lg.st KAt minutes. For example, in Figure 5-6, the alignment error at 25 
minutes is eight seconds of arc based on data taken between 20 and 25 mmutes. The phase 
of the low frequency motion is indicated with a sine-wave (dashed line). The open circles 
and triangles indicate data points from other phases shifted by a multiple of the basic 
period. 

The covariance matrix converges m two or three iterations to an asymptotic value. Hence 
the error does not converge. The uncertainty of the motion during the measurement mter- 
val KAt is dominant over the other noise inputs. 

The maximum error occurs when the rotational velocity is largest; the error is almost 
periodic with a period that is one-half the period of the motion. For At ^ 30 seconds the 
error is approximately proportional to jsin(ujt +<^q)| + ^. The residual error ^ is the 
consequence of other noise sources and quantization. These graphs illustrate that the rota- 
tional motion from the environment is the dominant source of error. For At> 30 seconds, 
the mstantaneous estimate is uniformly better than the simple average. On the other hand, 
for At^ 30 seconds, the instantaneous estimate has sporadic error spikes. Also for the 
At'^ 30 seconds, the effect of the LF motion is comparable to the effect of quantization and 
other noise sources. For the low frequency noise with a 29 mmute period, the results are 
graphed in Figure 5-9, 5-10, 5-11, and 5-12. Note that the error is not reduced by one- 
half even though the motion is one -half the previous value. These curves have the same 
characteristics as the 58 minute curves. 

The iterative alignment technique requires more computer computation and memory than 
the non-iterative techniques. In each application the increased computation requirements 
must be balanced with the gam in alignment accuracy to obtain an "optimum" alignment 
system. In some applications, the computer capability may prohibit using the iterative 
technique. The computation requirements can be reduced by using the assmaptotic value 
of the covariance matrix. 

5.3.1. 3 Gsnrocompass Alignment 

The alignment error for the simple average technique is graphed m Figure 5-13 as a 
function of alignment time, KAt. This error is the rms error based on ten trials. Since 
gyro quantization error is large, the alignment error is strongly dependent on the quanti- 
zation error in the east gyro output. In the simulation output we confirmed that al^nment 
error is directly correlated to the gyro quantization in the east direction. 
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If the instantaneous estimate of gravity is used, we will obtain an improved estimate of the 
gravity components, as indicated on Figure 5-13; compare graphs for K=1 and K=5. The 
Gyrocompass Alignment error will not be reduced as much as the level alignment error in 
Figure 5-13 since gyro quantization is dominant. 

5.3.2 Sensor Quantization 

Several tests were performed to determine the effect of sensor quantization. The results 
are in Table 5-3. The entries are the rms alignment errors based on ten trials. The 
following assumptions were made: the low-frequency rotational motion was one minute 
of arc with a one hour period; the first data input to the Estimation Routine was taken when 
the body frame was level. Orientation Case 1 was used; the word length was not restricted 
in the estimation routme. Also, K=5 for all tests and At = 30 seconds for one test and 5 
minutes for the other. A gaussian HF rotational noise was used. 

The results for level alignment are in columns A and C. Note that the posterior -mean 
estimate is sensitive to the quantization level for small values of At. Also, the accuracy 
is not improved when the quantization is reduced below 1.27 x 10“^ ft/sec. and 1.22 x 
10“^ rad. 

The results for Gyrocompass Alignment are in column D. The largest entries (116 for 
At =: 30 seconds and 120 for At = 5 minutes) result from a single trial in which there is 
one quantum error in the east-west gyro output. Note that even with the lowest quantiza- 
tion level (and At = 30 seconds) the quantization error in the east gyro is very significant. 

5.3.3 Computer Word Lergth 

As a first cut in determining the efforts of computer word length on alignment accuracy, 
a word length truncation was placed in the estimation routine in the computer simulations 
performed. A complete study of this effect could not be accomplished because of time 
limitation. The simulations were based on the same assumptions used for sensor quan- 
tization (Section 5.3.2) with quantization of 1.27 x 10 ft/sec. and 1.22 x 10 rad. 
Simulations were performed for word lengths of 15, 24, and 27 bits. The results of these 
simulations are shown in Table 5.4. The entries are the rms rotational alignment errors 
based on ten trials. Accuracy of the Level Alignment techniques are presented m columns 
A and C. The accuracy of the gyrocompass technique is presented in column D. 
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TABLE 5-3 

ALIGNMENT ACCURACY VS SENSOR QUANTIZATION 

Rotational Alignment Error (secj 

for • Gaussian Noise Distribution 
• Orientation Case I 
o K = 5 

A - Simple Average, Level 

C - Posterior- Mean Estimate of lastantaneous Components, Level 
D “ Simple Average, Gyrocompass 

At = 30 sec 


QUANTIZATION 

TECHNIQUE 

FT/SEC 

RADIANS 

A 

C 

D 

1.27 X 10"^ 

1.22 X 10"^ 

8,3 

3.7 

116.0 

1.27 X 10“^ 

1,22 X 10"^ 

8,4 

4.2 

10.2 

1.27 X 10'^ 

4.88 X lO”"^ 

9.0 

26,8 

10.5 


At = 5 min 


QUANTIZATION 

TECHNIQUE 

FT/SEC 

RADIANS 

A 

C 

D 

1.27 X 10"^ 

1.22 X 10"^ 

59.8 

33.3 

87.2 

1.27 X 10"^ 

1.22 X 10"^ 

60.0 

34.0 

120.0 

1.27 X 10'^ 

4.88 X 10"^ 

59.7 

35.5 

60,0 
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TABLE 5-4 

ALIGNMENT ACCURACY VS COMPUTER WORD LENGTH 
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The reader should be aware that the data given in Table 5. 4 is not conclusive for a deter- 
mination of the effects of word length on alignment accuracy. Only the estimation portion 
of alignment was considered. The program does not simulate the pre-processing or the 
determination of the T matrix. Further, the truncation performed was a simple truncation 
rather than a roundoff. Many of the computers that may be used with the ISU have roundoff 
capability and will be able to achieve better results than is indicated by the data in Table 5-4. 
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SECTION 6 


SECTION 6 

ALIGNMENT ACCURACY VERSUS CALIBRATION ACCURACY 


In the previous sections we considered the al^nment accuracy as a function of instrument 
noise, environment noise, instrument readout quantization, estimation technique and align- 
ment matrix computation. Referring back to Section 3, this leaves the consideration of 
the effect of calibration accuracy on the aligiunent accuracy. This invest^tion is the 
subject of this section. 

In Section 6. 1, the generalized error equations which relate calibration errors to errors 
in alignment are developed. In Section 6.2, these equations are evaluated for the orienta- 
tion cases, 1 to 4. In Section 6.3, worst-case alignment errors are developed from a 
worst-case combination of calibration errors for all four orientations. Finally in Section 
6. 4, the statistics of alignment errors are derived in terms of the statistics of the calibration 
errors. 

6. 1 GENERALIZED ERROR EQUATIONS 

The generalized error equations are developed in Chart 6-2 for both level and gyrocompass 
alignment. The assumptions leading to these equations are also listed in the chart. The 
Pre-Processing Computations, from which the error equations evolve, are presented in 
Chart 6-1. 

6.2 ERROR EQUATION FOR CLASS 1 TO 4 

The generalized error equations of Chart 6-2 are evaluated for each of the four orientation 
cases in this subsection. The results are listed in Charts 6-3 through 6-6. Note that the 
errors have been normalized by g and w and that A(Cjq) and A(Qjg) have been dropped 
since they have a negligible effect on the results. 

6.3 WORST-CASE ALIGNMENT ERRORS 

The results of Section 6.2 are evaluated in this subsection for the worst-case combination 
of calibration errors. These results are presented in Chart 6-7. The expressions in this 
chart are then substituted into the generalized error equations of Charts 4-5 and 4-6 to 
produce equivalent AT and ATT matrices, Charts 6-8 and 6-9. These latter results can 
then be evaluated, applying the results of Chart 4-3, Alignment Precision, to produce the 
basic measures of alignment errors, cone angles and total rotation angle. 
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CHART 6-1 


PREPROCESSING COMPUTATIONS 


Inputs ® 


The outputs Jj. 
computations 

Level GC 


{“•^)dt and (a*Bj^)dt (k = 1, 2, 3) are given by the following 


• - 

• "S|(EnT) 

• 

• Cfa.Aj^)^tl -pA/(Dj)j^-(Do)j^At 

. (^j^) = C(^.Gj^) Atj/At 

• fe’Aj^) s [(a-Aj^) AtV'^t 

~0 0 l1 nfaT^)" 

• (aT^ =10 0 {^) 

1 0 0 (aT^) 

“0 -1 0l RaT^" 

• (a^) =00-1 (aTSp 

0 10 &Ta^) 

, j.^t+At(^.g^)jjt ^ [ (u-Gj^)At }- [ (Bi)jj(a^)+(BQ)j^(iT^)+ (Bg)j^(^) JAt 

- L (Cii)k(^)MCsg)k (i^l^ 

*■ ^ ^ ‘"IS^k^-* -k^ ^-‘ -k^ ■*■ ^^OS^k^-’ 5k^ ^k^ ■*■ ^*'IO^k^-‘ -k^ ^-’ -k^ 

** C (Qji^k^ — * — k^ *** ^®IS^k^— ' — k^ (^*Sj^) - At 
vf^^‘(a«Aj^)dt = [(a-Aj^)Ati - (D 2 )i^(iT^)V - 0 ^)^^ (aT^)®At 

.J^t+At(^.^)^^ = j'^*’"^'(w.G,)dt 


1 "(G^-B2) 

Q^= -(Gg-B^) 1 -(Gg'Bg) 

-(G3. B^) -(Gg- Bg) 1 


• Q^= -(A,-B. 


■(Ai-Bg) -(A^-Bg) 
1 -(Aa-Bg) 


-(Ag-B^) -(Ag.Bg) 
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GENERALIZED ERROR EQUATIONS 


With the assumptions 

• All calibration constants are in error 

• The Q matrix errors are deviations from an Identity matrix. 

• No noise. 

• No quantization errors . 

P 

• S.- constant, = it = constant 

• As error coefficients Bj^) and (g. Aj^)=(g- Bj^,) 

These assumptions lead to the following error equations 

Level GC 

where 

• A(A4>)i/(A%(«^.B^) - ARj^ - A(Bj)^(£.B(^) 

- A(Bo)^(g.Oj^) - A(Bg)j^(g.Sj^) - A(Cii)i^(g. B^)2 

- ^(Cgs)^(g.Sj^)2 - A(Cjg)^(g.B^)(g.S^) 

- ^(Cog)j^(g.Oj^)(g. S^) - A(Cjo)j^(g. Bi^Hg- Oj^) 

. . A(g.A^) = -[ (g.B^) - (Dq)j^] A(D^y (Dj)j^ - A (Dq)^ 

-A(D2yg.B^)2-A(D3y£.Bj^)3 

0 -A(Aj^-B 2) -A(Aj*Bg) 

• • AQ^ = -A(A2'Bj^) 0 -AtAg'Bg) 

-A(Ag.B^) -ACAg.Bg) 0 

0 -A(G^ B2> -AtGj.Bg) 

• AQ*^ = -ACGg-Bj) 0 -A(G2*B3) 

-ACGg.Bj) -A(G3*B2) 0 

0 0 l1 r (g. B^) 

• (g-Ok)= 1 0 0 (S-§2^ 

1 0 oj [(g-B3)_ 

0-1 0l r (g. B^)' 

• (g-Sj^) =00-1 (g. Bg) 

0 1 oj j_(g Bg) 





CASE 1 ERROR EQUATIONS 


The nominal alignment is given by: 


[T] = 


0 0 1 
10 0 
0 10 


latitude = 45° 


The error equations are: 


Level GC 






E 


/ V <0® 


)-0.707(A(Gj-B2» -0. 707(A(Gj-B3)) 


' A(m= . BjA /a(A*)A /arA /A(Bg)2g\ ^(Cgg)2^^ 

U 0.7071 - -I 

/ V E / V 1 
/ \w / ^ w 




E 


KA®)2 


.E 


W 


E 


-0. 500 (a(Qjj) 2«^)-0. 707 (A(G^ • Bg)) 




.E 




^(A4>)g / ^ 


OJ 


E 


-0. 500 (A(Qjj)3to®) -0. 707 ( A(G^ . § 2 )) 


g / V g - 


|-(A(A^.:^)) 
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CHART 6-4 


CASE 2 ERROR EQUATIONS 


Th& nominal alignment is given by 


'0.577 0.577 0.577^ 

1= 0.815 -0.407 -0,407 
.. 0 0 707 -0.707 _ 


latitude = 45° 


The error equations are 


/a(A 4 >)i\ M(Bo)iS\ a 

a.40V^— j.0 57v\— 


(^^^sshA 


-0. 166(A(eijj)[j^)-0.907(A(Gi-B2))40,093(A(G|-Bj))- — 


/A(«A.B2)\ /A(A®)2\ /A(Bi)2B\ /^(Bo)2t'\ , 

U ^ j= 0 907l 1-0. 5771 — ~j-0 577 — ^*0.577' 


^^{Cn)2e^ h<=sshA h^ish^' 


-0.824(A(qjj)2'‘^-0.407(A(G2'Bj))-0.093(A(G2'B3))-- 


/A(A4.)g\ /'AlBjJje p(Bo)3M 

1- 0 093 ^0.577l /-0.5771 j^o.sn' 

‘ r \^E / V / 


/^(Egg)3gA A(Ejg)3!S^ 


1-0. 333^ ^ ^0.333^ ^ 


-0 00e(A(Qy)3ME)-O 407(A(G3 -Bj))-0.907(A(G3-B2))-- 


A(g-Bi) 


/A(D„) 


(Di)i/ \ K 


O.333(A(Dg)^tf)-0. 192(A(D3 )jB 


-0. 577(A (A j. Bg)) -0. 577 (A(a 1 - ^)) 


A(£.B2> 


^i(Dj)2\ MDoIz 


(Di)2/ V B 


0.333(A(D2)2B)-0 192(A ( 03)28 


-0.577(A(^ Bi))-0.577(A(^.B3)) 


Afg.Bg) 


^A(Djt)3\/A(Do)3 


-0.6771 )-0. 333(A(D2)g!;)-0 192(A(D3)3g' 

\(Dj) 3 / \ g / 


-0.577(A(^ Bj))-0.577(A(^-B2)) 
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CASE 3 ERROR EQUATIONS 


The nominal alignment is given by 


0 -0, 707 -0. 707 

[ T] = 1 0 0 latitude = 45° 

0 0.707 0.707 


The error equations are 


Level GC 






-0, 500\^ y-0.500^- — — 


L(A(g^.^)) 


0*2 / 


ACteE.^A A^\ 


0.7071 i-0.707 


^^(Cgs)2g^\ 


(A(G,-B,)) 


A(^^.Bgj\ /a(A$) 3 UaR 3 \ /A(Bj) 3 g\ ^( 33 ) 3 ^ 






A(g.Bj)\ A(Di)l\ A(Dq) 


Pi)i / \ g 


|+0.707(A(Aj^.B2))-0.707(A(Aj B 3)) 


0 607r^l f— 

g / ■ “•®%Dj) 2/A g , 


Lo.500(A(D2)2g)-0 354(Ap3)2gT 


-0. 707(A(.^.B3)) 


^(£•13) ] poPs] p°o)s 

g / ■■ ■ \(Dj) 3/ V g , 


1-0. 500(A(D2)3g)-0. 354(A(D3)3g2) 


+0. 707(A(.^.B2)) 
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CHART 6-6 


CASE 4 ERROR EQUATIONS 


1 ho nomuial aliyouncnt is {'ivon by 

0 407 -0.093 0 907 

I'T = 0 815 -0.407 -0 407 ialitudt - 45 

0 407 0 907 -0.093 

The crroi cqu4Uo»!> are 
Level GC 

\ A(a«)A /arA A(Bi)ib\ A(Bo)ik\ 


(At), > V' 


0 40?l -0 907\ -0 092| 


A<^05)jg' 

-0 166[ )-0 009 1-0 038( 1-0 Q84 

V O..E / V / \ / V „.E 


-0.333(A(Qb)u®)-0 577(A(Gj B^D-O 577(A(Ci B^)) 


I 0 577( 

. ! \(A«)2 ' 


/ARjX /A(B,),A A^sh«\ 

J 093 1-0 407 ItO 907l ) 

V(A«)2A,e/ V / V „E / \ / 

MC^)A A<^osk^\ 

f j-0 830(-^^0 034^ j-0.3B8l-_j 


-C 333 (A(Q|j) 2 W ^)-0 577{A(G2 Bi))-0 577(A(^ B 3 )) 

M(e;E 53)\ /A(A4)g\ /AR3\ 

. j, 907^-— )-0.40,^— j-0 093(— J 




I 

1 AVVrtJ 1 

f \ / 

V 4.E / 

l**U UOCH 1 

' / 


-0.577(A(Q^)3U.E).0 577(A(G3 Bj))-0.577(A(Gj. Bj)) 

/A(g-Bi)\ /A(Di)j\ /A(D3 )j\ 

I 1 =-0 5071 ^^-1 1-0 166{A(D2)ig)-0.068{A(D3)jg‘^; 

' B ' Vi)/ ^ B ^ 

+0 033(A(Aj-B2))-0.907(A(Aj 53)) 

/A(s.Ba)\ /A(Di) 2 ' /A(Do) 2\ 

I 1 = .0.007( U -0.009(A(D2)2g) 

g ' 8 ' 

-0.407(A(^ Bj))-0.907(A(^ ^)) 


l3)\ 


1 8 ) 

' Ul>3 ) V 8 ) 


-0.407(A(^.Bj))-0 407(A{^ Bj)) 
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WORST CASE CALIBRATION ERRORS 



Case 1 

Case 2 

Case 3 

Case 4 

2E^ + I.4IE2 

4.06E^ + 1.57Eg 

3a41Ej^ + Eg 

2,70Ej + 2 .O 7 E 2 

3Ej^ + 1.' 66Eg 

4. 06Ej + 2. 23 E 2 

3« 91Ej|^ + ^2 

3'69Ej^ + 2 .O 7 E 2 

+ 1, 66E2 

4, 06E^ + 1. 42 E 2 

3a 91Ej^ + 2 E 2 

3« 36E^ + 2o O6E2 

1* 

2. 36Eg + E^ 

1.51Eg + E^ 

1.74Eg + E^ 

1, 

2a 36Eg + E^ 

2. 17Eg + E^ 

1« 52E2 + 

3o 

2. 36Eg + E 4 

2» STEg 

3 . 39Eg + E^ 


Gyro errors: 


Accelerometer errors. Eg = 


AR ABg I ACg'^ 





A(A. B) = ADgi 



, E, = 


B) 

= 

AQ w® 


ADq 



g 
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WORST CASE LEVEL ALIGNMENT FROM CALIBRATION ERRORS 







Case 1 

Case 2 

Case 3 

Case 4 



ATji 

1* 1 Eg 

3« 14 Eg +1.33 E^ 

1,51 Eg +E 4 

2.77Eg +1, 24E^ 

ATj2 

1. 1 Eg + E^ 

3 . 14 Eg +1.33E^ 

2. 27 Eg +E 4 

1.85Eg +1. IIE 4 


0 

3, 14 Eg +I. 33 E 4 

2. 27 Eg +E 4 

1.37 Eg +0. 63 E 4 

^^21 

0 

2. 22 Eg +0. 94 E^ 

0 

1. lOEg +0. 5 OE 4 


0 

5. 55Eg + 2 . 36 E 4 

1.07 Eg +O. 7 IE 4 

1.07 Eg +0. 59 E 4 

^"^23 

1. lEg +E 4 

5. 55 Eg + 2 . 36 E 4 

1.07 Eg +O. 7 IE 4 

3. 19Eg +1. 53 E 4 


0 

0 

0 

0. 55Eg +0. 25 E 4 

^^^32 

0 

5. 78 Eg +2.45E^ 

2. 27 Eg +E^ 

0. 3lEg +0. I 6 E 4 

^’^33 

1. 1 Eg + E^ 

5.78Eg +2.45E^ 

2. 27 Eg +E^ 

2. 50Bg +I. 44 E 4 



(ATTT)j,2 

!• 1 Eg + E^ 

3. 85 Eg +1. 63 E^ 

1. 51 Eg +E 4 

2. 84 Eg +1. 63 E 4 


(ATTT)^3 

1. lEg + E^ 

1 

7.70Eg +3.27E^ 

2. 14 Eg +1.4. E^ 

2.13Eg +I. 4 IE 4 

(ATTT>23 

0 

3. 14 Eg +I. 33 E 4 

0 

1.66Eg t 0 . 49 E 4 
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WORST CASE GYROCOMPASS FROM CALIBRATION ERRORS 



Case 1 

Case 2 

Case 3 

Case 4 

AT,, 

1 * lEg + 

3.13Eg + I. 33 E 4 

I. 5 IE 3 + E 4 

2o73Eg + 1 « 24E^ 


1 , 

3.13Eg + 1.33E^ 

2 « 27Eg + E^ 

I. 85 E 3 + l.llE^ 


0 

3, 13 E 3 + 1. 33 E 4 

2 . 27Eg + E^ 

1 . 37 E 3 + 0 . 63 E 4 

ATg, 

0 

2 . 22 Eg + 0 . 94 E 4 

0 

2. 94E, + 1. 95Eg 

^^22 

1.41E, + Eg 

+ 1 * lEg -i- E^ 

6 . 64E, + 2 . 77Eg 
+ O. 5 OE 3 + 0.21E^ 

3.41E, + Eg 
+ 2.14Eg + 1.41E^ 

6 . 53E, + 4. 33Eg 
+ 3 . 42 E 3 + 1 „ 63E^ 

^^23 

1 . lEg + E^ 

6.64E, + 2.77Eg 
+ 4* 93Eg + 2» lOE^ 

3.41E, + Eg 

0« 66Ej^ 4“ 0» 44Eg 
4 * 3* 42Eg + 1. 63E^ 


1.41E, + Eg 
4 - 1. lEg + E^ 

7 . 65E, + 3. 20Eg 
+ 3.49Eg + 1.48E^ 

4.83E, + 1.41Eg 
+ Id 51Eg + E^ 

5 . 87E, + 3. 89Eg 
+ 2 . 77 E 3 + 1.24E^ 

^^32 

0 

3.83E, + 1. 60Eg 
+ 2 . 7 IE 3 + I.I 5 E 4 

6 . 8 OE 3 + 3E^ 

2. 94E, + 1.95Eg 
+ 1 . 37 E 3 + 0 . 63 E 4 

^'^33 

1. lEg + E^ 

3.83E, + 1. 60Eg 
+ 2 . 7 IE 3 + I.ISE^ 

2. 27Eg + E^ 

2.94E, + 1.95Eg 
4- 1, 85Eg 1. IIE^ 


(ATT’’^),g 


3.84Eg + 1.63E^ 

I. 5 IE 3 + E^ 

3 . 42 E 3 + 1.63E^ 

( ATt '^),3 

l.lEg + 

3. SSEg + I. 4 IE 4 

3 . 2 IE 3 + I. 4 IE 4 

2 . 4 OE 3 + I. 4 IE 4 


I 9 41 E^ + E^ 

4 - 1, 1^^ 4 - E^ 

9. 38Ej 4 - 3. 92 E 2 
+ 3* 84Eg 4 - 1. 63E^ 

4.83E, + 1.41Eg 
4 " l*51Eg 4 " 

7 . 20E, + 4.78Eg 
4 * 3»42Eg 4 - 1. 63E^ 
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6.4 STATISTICAL ERRORS 

A set of one-sigma values of the calibration coefficients was provided, by ERG and are 
listed on Chart 6-10. These values were used to derive the one-sigma values of 


o B ) 
-k 


by applyn^ the results of Charts 6-3 throi^h 6-6 of Section 6,2. The one-sigma values of 
these quantities, defined as cr^j^ andcr^j^ respectively, are listed in Chart 6-10 for all four 
orientation cases. Finally, these latter values along with the results of Section 4, 3 were 
used to derive the one -sigma values of the elements of the At and ATT matrices for both 
level and gyrocompass alignment and for all four orientations. These values are presented 
in Charts 6-11 and 6-12. 

An analysis of the statistical calibration errors listed in Chart 6-10 showed a strong 
dependence of the results on the particular one-sigma value selected for gyro bias. Addi- 
tional statistical analysis of this coefficient is indicated. The development in this section 
has been presented in enough detail to permit a calculation by the reader of statistical 
alignment errors for a different set of coefficient one-sigma values, should this be desirable 
in the future. 
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STATISTICAL CALIBRATION ERRORS 



Case 1 

1 

Case 2 

Case 3 

Case 4 

*^wl 

12.0 X 10”® 

12. 1 X 10”® 

12.2 X 10”® 

12.0 X 10"® 

‘^w2 

12-3 X 10”® 

12.1 X 10”® 

12.2 X 10"® 

12.3 X 10"® 

^uiZ 

12.3 X 10”® 

12. 1 X 10”'® . 

12.2 X 10”® 

12.2 X 10”® 

'^gl 

12.1 X 10”® 

12.5 X 10”® 

12.1 X 10”® 

12.4 X 10"® 

<^g2 

12.1 X 10"® 

12.5 X 10"® 

11.5 X 10"® 

11.8 X 10"® 

>^g3 

12.9 X 10”® 

12,5 X 10”® 

12.7 X 10”® 

13.4 X 10"® 


Gyro Bias 

Scale Factor 
Unbalance 
Compliance 
Misalignment 


Coefficient One-Sigma Values 


0. 1 deg/hr 
10 "'^ 

0. 15 deg/hr/g 
0. 04 deg/lir/ 

3 X 10“^ 


Accelerometer Bias 
Scale Factor 
Second Order 
Third Order 
Misalignment 


6.7 X 10"® g 
10 "® 

0. 9 X 10”® g/g^ 
0. 1 X 10”® g/g^ 
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STATISTICAL LEVEL ALIGNMENT FROM CALIBRATION ERRORS 



Case 1 

Case 2 

Case 3 

Case 4 

‘^ATll 

12. lx 10"® 

10.2 X 10"® 

12. lx 10"® 

11. 5 X 10"® 

*^AT12 

12. lx 10"® 

10.2 X 10"® 

8. 6 X 10"® 

11. 8 X 10"® 

“^ATIS 

0 

10.2 X 10"® 

8. 6 X 10"® 

5.3x 10"® 

°’at21 

0 

7.2 X 10"^ 

0 

4. 6 X 10"® 

^^AT22 

0 

9.5 X 10“® 

8. 6 X 10"® 

5.4x 10"® 

‘^AT23 

12. lx 10“® 

9.5 X 10"® 

8.6x 10"® 

11.5 X 10"® 

"^ATSl 

0 

0 

0 

2.3 X 10"® 

^^At32 

0 

10. 8 X 10"® 

8. 6 x 10"® 

1.3x 10"® 

‘^AT33 

12. 1 X 10"® 

10. 8x 10"® 

8. 6 X 10"® 

12.8 X 10"® 


^ATTT12 

12. lx 10"® 

12. 5 X 10"® 

12. 1 X 10"® 

12.5 X 10"® 

°"ATTT13 

12. lx 10"® 

12. 5 X 10"® 

12. 1 X 10"® 

11.9 X 10"® 

^ATTT23 

0 

8. 8 X 10"® 

0 

5.3 X 10"® 
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STATISTICAL GYROCOMPASS FROM CALIBRATION ERRORS 



Case 1 

Case 2 

Case 3 

Case 4 

CT^Tll 

12.1 X 10“® 

10.2 X 10"® 

12.1 X 10"® 

11.5 X 10"® 

^AT12 

12. 1 X 10“® 

10.2 X 10"® 

8.6x 10"® 

11.8 x 10"® 

<^AT13 

0 

10.2 X 10"® 

8.6 X 10"® 

5.3 X 10"® 

°'aT21 

0 

7.2 X 10"® 

0 

7. 0 X 10"® 

o'AT22 

17. 0 X 10"^ 

12.1 X 10"® 

12.2 X 10"® 

15.5 X 10"® 





°'AT23 

12. 1 X 10"® 

12.1 X 10"® 

12.2 X 10"® 

1.6 X 10"® 




°"aT31 

17.0 X 10"® 

14. 0 X 10"® 

17.3 X 10"® 

14. 0 X 10"® 




' 

<^AT32 

0 

7.0 X 10"® 

8.6x 10"® 

7.0x 10“® 



"^ATSS 

12.1 X 10"® 

7. Ox 10"® 

8. 6x 10"® 

7. 0 X 10"® 




^ATTT12 

12.1 X 10"® 

12.5 X 10"® 

12.1 X 10"® 

12.5 X 10"® 

^ATTT13 

12. 1 X 10"® 

12.5 X 10"® 

12.1 X 10"® 

11,9 X 10"® 

a 

ATTT23 

17. Ox 10"® 

17.1 X 10"® 

17.3 X 10"® 

17.1 X 10"® 
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APPENDIX 



APPENDIX A 

SUPPORTING CALIBRATION NOISE ANALYSIS 


Section 2<,2 develops the equations relating coefficient errors to sensor output error, AP, 
for both gyros and accelerometers. Charts 2-3 and 2-4 list these equations together with 
equations which describe ^ terms of the statistics of the noise environment. 

The curves in Figures 2-1 through 2-4 show how these coefficient errors vary with the 
tune interval. At, over which the sensor output is averaged. 

In this Appendix we show the development of the noise equations ui more detail than was 
given in Section 2.2. Section A. 1 outlines the approach to the problem and relates the 
statistics of the output to those of the input to the process. The stochastic process con- 
sidered includes the geometry of the problem and the averagu^ as well as the sensor it- 
self, for the AP of Section 2.2 is the sensor output after it has been averaged over the 
time interval. At. In fact, sensor characteristics play a relatively minor role in the 
noise analysis. 

Sections A. 2 and A. 3 discuss the geometrical effects and the statistics of the sensor mput 
for gyro and accelerometer respectively. Finally, Section A. 4 gives some of the mathe- 
matical details and approximations employed in the calculations for Figures 2-1 through 
2-6 in Section 2. 2. 

A. 1 APPROACH TO THE PROBLEM 

The data, P, used in calibration is the sensor output summed over a time interval, At. 
The noise-induced error is Ap and it is our aim here to relate AP to the input noise. 

The sensor, together with the subsequent summation of the pulses, is considered a 
stochastic process. The input to this process is a stochastic variable (the environment 
noise) plus, m the case of the gyro, an mput equivalent to the internally-generated noise. 
The process model adopted for both gyros and accelerometers is 



AP 


A-1 





where An is the input equivalent sensor noise (gyro only), and Aw*G and Aa.A represent 
environment noise inputs to gyro and accelerometer, respectively. The summation has 
been replaced with an integration which means we are neglecting the sampling process and 
replacing it with a continuous process. This has little effect' on the results since the 
sample interval is small compared with At. 

Then the power spectral densities* of input and output are related by: 


Pp(f) = |Ts(iu>)|2|Tj(jco)|2p^(f) 


(A-1) 


where Pp(f) 

Tg(jw) 

Ti(]to) 


power spectral density of oulput error, Ap 
power spectral density of the input 
transfer function of the sensor 
transfer function of the integration. 


The variance of the output error is related to the output power spectral density by 
(assuming stationary noise): 


var [AP] = J* P ff)df (A -2) 

-CO P 

Finally, the rms error is the square root of the variance. 

The following sections discuss the gyro and accelerometer separately. As will be shown, 
the effect of environment vibrations upon both types of sensor depends on the sensor's 
orientation wrt** environment noise and also wrt^ and g, respectively. For this rea- 
son two sensor orientations are investigated (vertical and horizontal) and simplifying 
assumptions are made concerning the orientation of environment noise wrt these directions. 

It is appropriate to discuss the integration transfer function Tj(]w) at this pomt smce it 
is common to both types of sensor. 


♦Other approaches might be taken to the problem. This approach is dictated by the 
nature of the available noise data. 

**wrt = with respect to. 
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Let AP(t) be 62 g)ressed as 


m 

AP = J x(T)w(r)dT 

- 


where x(t) is the iaput to the integration and w(t) is a boxcar weighting function defined 
as one in the interval t-At ^ t s t and zero outside this interval. Then the transfer 
function is: 


At 


Tj(]w) = LCw(t)] = J e 


-jwt 


dt 


,-jwAt 


-1 




j (li = 2iri 


and the square is; 



2 


(gjcoAt_ij 

(-]w) Oco) 


2 (1-cos coAt) 



A. 2 GYRO 

The input sensed by the gyro is w . G 

where G = a unit vector along the gyro input axis 

w = the total angular velocity of the gyro with respect to inertial space. 

Now the gyro is on a turntable rotating with a constant angular velocity with respect 
to the laboratory. The laboratory (the "environment") is subject to vibrations due to 
traffic, earth tremors, 'wind, etc. This is represented by saying the laboratory oscil- 
lates with an angular velocity w® with respect to inertial space. Hence 

T e E 
w = w +w +w 

and w*G = w'^»G + G + G 
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Then the error in w . G is 

A(ti)*G) = A(oj'^*G) + A(co®»G) + A(co^»G) 

This can be simplified by observing that the first term on the right is zero if the gyro is 
firmly attached to the turntable. 

Let w® = w® + Aw® in the second term where w® is the mean of w®. The mean has been 
assumed zero, so that we need only be concerned with Aw®, the distribution about the 
mean. Then the second term becomes Aw®*G. To simplify the problem and to get a 
''maximum rms" error, we will neglect the wT modulation of Aw®-G and compute the 
errors due to environment angular displacement noise for the two cases: gyro vertical 
and gyro horizontal. 

For the A(w‘'^»G) term we have to consider in what coordinate system these omegas are 
defined. The dot product is invariant with respect to a bnear coordinate transformation, 
but to be meaningful both vectors must be expressed in the same coordinate system. For 
our purposes here it is a matterof indifference whether G is fixed or w is fixed. We 
want to focus attention on the angle between them, which varies due to the vibration of the 
laboratory environment. (The angular velocity noise Aw® is the time derivative of this 
angle.) Let 0 =0^-f A0 be this angle. Then 

w^*G = w^cos 9 

and A(w^*G) = w^A cos 0 

= {-w®sin0Q)A0 

where A0 can be termed the angular displacement noise. Thus an important part of the 
input noise error arises from the gyro sensing a variable component of earth rate, the 
variation being due to vibration of the surrounding environment. 

Putting all together we find that the mput error is 

A(w*G) = Aw®*G + (-w^sin0Q)A0 
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Finally, the input power density spectrum is 


Pi„W - 


input spectrum equivalent ■) 
of sensor noise j 


Power spectrum 1 
of A{w*G) ) 


" ^n^^) + cos^0^P^e(f) + (w®sin0Q)^Pg(f) 


(A-3) 


where P (f) 

IV ' 

P„e(f) 

Pe(f) 


power spectral density of gyro internal noise 
power spectral density of angular velocity noise Aw® 
power spectral density of angular displacement noise, A0 
angle between G and the horizontal plane. 


The angle 0^ enters via the dot product Aw®*G. For completeness, there should be an 
azimuth angle here too, but this dot product is independent of azimuth since we assume 
the angular velocity noise is isotropic in the horizontal plane and zero in the vertical. If 
it were completely isotropic we could drop the cos0^ factor also. 


The variance of the gyro output noise is found by substituting (A-1) and (A-3) into (A-2) 
and performing the mdicated integration, after functional expressions are found for the 
transfer functions and power spectra. It turns out the resulting integrand is not integrable 
in general, so the integration is performed numerically. 


The angular displacement noise spectrum Pg(f) is that given in Figure 3-4, Section 3.2.2, 
of the Development Document. This must be converted to (radians)^ rms/cps before being 
used in the computations. The angular velocity noise spectrum, P^e(f)j is derived from 
the same figure by multiplying by (2irf)^ = (juj)(-]w). This corresponds to differentiation 
in the time domain. Since the angular displacement (and hence the angular velocity) is 
assumed isotropic in the horizontal plane and zero in the vertical, cos0q m Equation (A-3) 
is equal to zero for a vertical position of the gyro-sensitive axis and equal to one for a 
horizontal position. Furthermore, sinfl^ = 1//T for both orientations, since 0^ is 
assumed equal to 45*^. 


To summarize, we have for the variance of the gyro noise error 

00 

var CAP] = J P„(f)df 

—CO P 
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where Pp(f) = 1 Tg(jw) | ^ | Tj(]w) 1 ^ 




P^{f) A 


Pfl(f) 


, for gyro vertical 


= |Tg(iw)|2|Tj(]w)| 


'OJ 




Pn(f) + — 1 P0(f)+Po;eW 


, for gyro horizontal 


Sensor dynamics for the gyro was approximated hy a 20 Hz first order loop. Consequently, 


Ts(jw)i 


I 2 • 


1 + (f/20)" 


is adopted for the squared transfer function. 


A. 3 ACCELEROMETER 

The description of accelerometer input is similar to that of the gyro. The input sensed 
by the accelerometer is a • A , where 

A = a unit vector along the input axis of the accelerometer 
a = total acceleration 
= a® +g 

a® = environment acceleration (noise) 
g = acceleration of gravity 

The unwanted input, or noise, is 


A(a • A) = A(a® • A) + A(g . A) 


Assume: Ag 


0 

+ Aa® , with the mean, a® , equal to zero. 


Then 


A(a • A) = Aa® • A + A(g • A) 



The first term on the right Is the component of environment noise along A and the second 
term is the variation in the component of g sensed by the accelerometer; the variation 
being due to the angular vibration of A with respect to g. H 0 = 6 ^ + A 8 is the angle 
between A and g, then: 


A(g.A) = (-gsin 0 Q)A 6 


Then the input power spectral density is 

P^(f) = cos 20 ^p^e(f) + (gsin 9o)^Pg{f) (A-4) 


where Pg^eCf) 


S 


power spectral density of the environment acceleration noise 
power spectral density of the angular displacement noise, A 0 
angle between a® and A 
nominal angle between g and A 
magnitude of g 


The variance of the accelerometer output noise is found by substituting (A-1) and (A-4) 
into (A-2) and integrating numerically. 

The noise spectra P^^eW ^ 0 (^) those given in Figure 3-3 and 3-4 in Section 3,2.2 
of the Development Document. Figure 3-4 gives Pg(f) in sec.^ rms per cps. 

Since the random acceleration inputs are assumed isotropic, cos0^ is equal to 1 in 
equation A-4, for all orientations of the accelerometer. But since the random rotational 
inputs are assumed isotropic only in the horizontal plane (and zero about a vertical axis), 
sin0Q in equation A-4 is one for a horizontal accelerometer and zero for a vertical 
accelerometer. 


To summarize then, we have 

var CAP] = J“p (f)df 

-CO P 

and Pp(f) = |Tg(]iu)| ^|Tj(]w) ^P^e(i) ; for acc. vertical 

~ I Tg(]w)l ^ |Ti(]w) ^ [Pae(f) + g^Pg(f)] , for acc. horizontal 
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Sensor dynamics for the accelerometer was approximated by a 1000 Hz first order loop. 
Then the squared transfer function is 

IXgaco)! 2 = ! 

1 + (f/iooor 

A. 4 COMPUTATIONAL DETAILS 


This subsection gives the computational details and approximations used in the computa- 
tions for the noise curves in Section 2o2, Figures 2-1 through 2-6. Some of the approxi- 
mations are rather rough, but this is not unreasonable in view of the approximate 
character of the power spectra and the fact that these spectra may not be representative 
of the actual noise environment at the Laboratory site. The results should be interpreted 
with suitable reservations. We first discuss the gyros, then the accelerometers. 

Gyro 

For the gyro we have 

CO 

var [AP] = J |Tg( 30 J)p|Tj(jto) 12 [pjf) + cos^ 0^P^e(f) + (‘^^sins/ Pg(f)3df (A-5) 

where |Tg(]w)P = 

1 + (f/20)^ 

, 19 2(1 - cos2nf At) 

; — 


and where cos 0^ = 1 for the gyro in a horizontal position 

= 0 for the gyro m a vertical position 
and sia0Q = 1//T in both positions. 

For the input equivalent to the internal noise, we adopted a spectrum which is flat up to 

-1 

a frequency of (100 hours) and which yields an rms noise of 5 x 10 degree/hour when 
integrated from (14 hours) to «>, and havii^ a peak at 16 cps. For computation, this 
spectrum was approximated piecewise as given in the following table: 
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Table A-1 


f(cps) 

((deg/hr)^/cps) 

0 to (100 hrs)"^ 

3. 02 X 10'^ 

(100 hrs)"^to 10"^ 

2.33 X 10"^^f^ 

10"^ to 8 

2.33 X lO"*^ 

8 to 24 

lx 10"® 

Above 24 

6 x 10 "^f"^ 


For the ai^ular displacement noise spectrum, Po(f), Figure 3-4 of subsection. 3.2.2 of 
the Development Document was used. This figure gives the spectrum in (deg/hr) per cps. 
When this is converted to (radians)^ per cps and multiplied by (w^/^/2)^ with = 

15 deg/hr, it is approximated piecewise as m the following table; 


Table A-2 


f(cps) 

( CO® sine^)^Pg(f) ((deg/hr)^/cps) 

0 to 5 X 10"® 

0. 106 

5 X 10"® to 0.317 

2.64x lO'^f"^ 

0.317 to 3.17 

2. 64 X 10"^ 

3.17 to 15.7 

2.64X 10"^® 

Above 15. 7 

1. 15 X 10"®f 


The angular velocity spectrum, which is needed only for the horizontal position, was 
derived from Table A-2 by multiplication by 

(]w)(-jw) = (27Tf)2 

which corresponds to differentiation in the tune domain. 

The integrand in equation (A-2) is not integrable due to the j Tj.(j(n)j^ factor, so this 
factor was also approximated piecewise. This factor is one at f=0 and descends to zero 
at f = 1/ At with a shape similar to cos (2 it f At). Beyond this frequency it proceeds in 
arcs of rapidly diminishing amplitude with zeros at 1/kAt, k = 3, 5, 7, .. . . For the 
numerical work, this was approximated by a constant 0. 7 from 0 to 1/At and by the 
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envelope of the peaks, 1/(17 fAt)^ , thereafter. Since the breakpoint depends on the 
calibration time. At, a new approximation is made for each value of At plotted. For 
the curves of Section 2.2, the noise errors were computed for calibration times of 
2, 5, 10, 15, 20, 30, 40, and 50 minutes. 

As mentioned earlier, the sensor transfer function |Tg(ja))]^ does not change the 
results significantly. 

It may be of interest to compare the separate contributions of the noise sources. In 
terms of contributions to the variance, and using the above approximation, we have, 
m (deg/hr)^ : 


„ 2.6x10“'^ 5.5x10"^ 

Contribution of sensor noise = 1. 7 x 10" ‘ + + 

At (At)2 

7 8.5x10"^^ 

Contribution cf A8 noise = 1.8x10" -2.4x10""^ At - 

(At)2 

4.35x10“^ 2.11x10"^ 

Contribution of w® noise = + 

At (At)^ 

Accelerometer 


For the accelerometer, we have: 


var C Ap] 



Ti(]w) 




sm2 0^Pg(f)]df 


(A-6) 


where lTg(j^)j^ = 

1 + (f/1000)^ 

{2TTif 

2 

as given previously. The factor g multiplying the Pg(f) term has been omitted because 
we want the variance in units of g^. Then the standard deviation will be in units of g. 

As with the gyro, Pg(f) from Figure 3-4 of Section 3.2.2 (Development Document) must 
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2 

be converted to (radians) /cps. The environment acceleration noise spectrum P„-(f) 

n <10 

from Figure 3-3 cf the same section is already in units of g /cps. 

The sensor transfer function is effectively set equal to one for the same reasons discussed 
above for the gyro case. The numerical approximations for the angular displacement 
spectrum, PqCI), and for the integration transfer function, jTj(ju)|^, have also been 
given in the discussion of the gyro. 

The environment acceleration spectrum, is from Figure 3-3 of the Development 

Document. For numerical work, this was approximated stepwise as follows: 


Table A-3 


f (cps) 

Pae(^) (g^rms/cps) 

0 to 10"^ 

JO-15 

10"® to 10"^ 

10"^^ 

10"^ to 1 

10"® 

1 to 10® 

10"® 

10^ to 10® 

10"^^ 

Above 10® 

0 


For the accelerometer in a horizontal position, the major source of error is the angular 

displacement noise which causes the accelerometer to sense a varying component of the 

2 

gravity vector. The separate contributions to the variance of AP are (in units of g ): 


Contribution of environment? ^ 
acceleration ) 


1.1 X 10"^^ 
+ 

At 


1.92 X 10"® 
(At)^ 


Contribution of g-pickup ? i. q ^ ^ in-9 

due to A0 noise » “ ^ 


2.3 X 10"^'^ At 


6.7 X 10"^® 
(At)2 


If a bubble is used to correct for the low frequency portion of the angular displacement 

noise, either mechanically or mathematically, the contribution of this error source 
2 

is (in units of g ): 
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Contribution of A0 noise ( 
based on Model 1 f 

= 7. lx 10"^^ - 

2. 3x10"^^ At + 

1 

Contribution of A0 noise 1 

1.5x 10"^ 

2.4x lO""^ 

based on Model 2 j 

At 

<1 


6» 7 X 10' 


(At)^ 


•15 


for the two modifications of angular displacement noise spectrum that were investigated. 
These two cases represent two models for the angular displacement spectrum after 
correction by the bubble level. In Model 1 the spectrum is simply assumed to be zero 
for frequencies below (50 min)”^. In Model 2 the spectrum is assumed to be: 

K 

= 

In this case the contribution to the variance of AP is: 


var 


At 

[J 0(t)dt] 

0 


2C, 


-CaAt 


,2 L 


■ 

1 + C2At 


where = a^^KTr/2 
For the computations we took 

2 2 
Cj = (4. 5 sec) (converted to radians ) 

C2 = 2it X 10“Vsec 
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